Effects of different geometries on the conductance, shot noise and tunnel 
magnetoresistance of double quantum dots 
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The spin-polarized transport through a coherent strongly coupled double quantum dot (DQD) 
system is analyzed theoretically in the sequential and cotunneling regimes. Using the real-time 
diagrammatic technique, we analyze the current, differential conductance, shot noise and tunnel 
magnetoresistance (TMR) as a function of both the bias and gate voltages for double quantum 
dots coupled in series, in parallel as well as for T-shaped systems. For DQDs coupled in series, we 
find a strong dependence of the TMR on the number of electrons occupying the double dot, and 
super-Poissonian shot noise in the Coulomb blockade regime. In addition, for asymmetric DQDs, we 
analyze transport in the Pauli spin blockade regime and explain the existence of the leakage current 
in terms of cotunneling and spin-flip cotunneling-assisted sequential tunneling. For DQDs coupled 
in parallel, we show that the transport characteristics in the weak coupling regime are qualitatively 
similar to those of DQDs coupled in series. On the other hand, in the case of T-shaped quantum dots 
we predict a large super-Poissonian shot noise and TMR enhanced above the Julliere value due to 
increased occupation of the decoupled quantum dot. We also discuss the possibility of determining 
the geometry of the double dot from transport characteristics. Furthermore, where possible, we 
compare our results with existing experimental data on nonmagnetic systems and find qualitative 
agreement. 



PACS numbers: 72.25.Mk, 73.63.Kv, 85.75.-d, 73.23.Hk 



T3 

C 

O 

i o i 

(N 
> 

m 
o 

^6 
o 

OO 

o 



X 



I. INTRODUCTION 

Transport properties of double quantum 
dots (DQDs) have recently attracted much 
interest . 1 - 2 - 3 - 4 - 5 - 6 - 7 - 8 - 9 - 10 - 11 ' 12 - 13 - 14 - 15 ' 16 - 17 - 18 - 19 - 20 - 21 - 22 - 23 - 24 

Since the behavior of double quantum dots resembles 
the behavior of molecules, DQDs are frequently called 
artificial molecules, and are thus considered as ideal 
systems to study the fundamental many-body inter- 
actions between single electrons and spins . 25 i 26 i 27 i 28 i 29 
Double quantum dots exhibit a variety of interesting 
effects, such as for example current rectification due to 
the Pauli spin blockade j^£ i 15 ' 16 i 17 ' 18 negative differential 
conductance' 



pumping fHiii Kondo effect ^-^L^l 



formation of molecular statesj^i 3 .^ S p m 
t) 30 1 3i 1 32 etc Furthermore, 
double quantum dots are also interesting for future ap- 
plications in quantum computingi 33 i 34 i 35 i 36 In addition, 
when the leads are ferromagnetic, transport properties 
of the system strongly depend on the relative orientation 
of the magnetizations of electrodes, leading to the tunnel 
magnetoresistance (TMR) effect, spin accumulation, 
exchange field, eic ^ML^^L^. 

The problem of spin-polarized transport properties has 
been so far mainly addressed in the case of single quan- 
tum dot; . 39 - 40 - 41 - 42 - 43 - 44 - 45 - 46,47 - 48 - 49 - 50,51 - 52 - 53 - 54 In partic- 
ular, in the strong coupling regime it was shown that the 
Kondo peak becomes split when the magnetic configu- 
ration changes from the antiparallel to parallel one, and 
that this splitting can be compensated upon applying ex- 
ternal magnetic fiekh 51 ' 52 i 53 In the weak coupling regime, 
on the other hand, the parity effect of the linear response 



TMR was predicted^ 3 - and the zero-bias anomaly was 
found in the differential conductance when magnetic mo- 
ments of the leads form an antiparallel configuration. 44 
In this paper we extend the existing theoretical studies 
and consider transport through coherent double quantum 
dots weakly coupled to external ferromagnetic leads. We 
note that spin-dependent transport properties of DQDs 
have been analyzed very recently in a few theoretical 
papers i 16 ! 23 ' 55 The considerations were however limited 
to either first-order tunnelin g or the Coulomb blockade 
regime. In particular, in Ref. I16| the formation of a pure 
triplet state was predicted, Ref. [23| deals with transport 
in the deep Coulomb blockade regime, while Ref. [HH 
presents a detailed analysis of sequential transport in the 
case of non-collinearly polarized leads. 

The goal of this paper is thus to analyze the transport 
properties of double quantum dots coupled to collinear 
ferromagnetic leads in the full weak coupling regime, 
i.e. including sequential, cotunneling and cotunneling- 
assisted sequential processes. Furthermore, we also ana- 
lyze the effect of different geometries of the double dots 
on transport, in particular, the cases of DQDs coupled in 
series, in parallel, as well as T-shaped double quantum 
dots are considered. The comparison of numerical results 
obtained for different geometries may in principle help in 
determining the system's geometry, which may be of im- 
portance in discussing and understanding experimental 
results, especially on self-assembled quantum dots. The 
present analysis is based on the real-time diagrammatic 
technique which consists in a perturbation expansion of 
the density matrix of the system under consideration, 



and the relevant operators, with respect to the coupling 
to external leads. The advantage of using the real-time 
diagrammatic technique is that it takes into account the 
effects of the exchange field and the renormalization of 
the dot levels in a fully systematic way. In particular, 
in the case of asymmetric DQDs coupled in series, we 
consider transport in the Pauli spin blockade regime and 
show that the leakage current in the blockade region re- 
sults from cotunneling and cotunneling-assisted sequen- 
tial tunneling processes. As far as the shot noise is con- 
cerned, we show that the noise is super-Poissonian in the 
Coulomb blockade regime and drops to sub-Poissonian 
value in the sequential tunneling regime, where the Fano 
factor approaches one half. This behavior is observed 
for parallel and serial DQD's geometries. In turn, for 
T-shaped double quantum dots, we find a large super- 
Poissonian shot noise due to increased occupation of the 
decoupled dot. We also analyze the TMR effect and find 
its strong dependence on the transport regime and num- 
ber of electrons in the DQD in the ground state. For 
DQDs coupled in series and in parallel the TMR takes 
the values ranging from around a half of the Julliere TMR 
to its full valuer 7 , whereas for T-shaped DQDs the TMR 
may be enhanced above the Julliere value. 

The systems considered in this paper may be realized 
experimentally for example in lateral and vertical semi- 
conductor quantum dot s 6 ' 8 ' 10 ! 56 ' 57 or single wall metal- 
lic carbon nanotubes with top gate electrodes i 12 i 58 ' 59 i 60 
The latter systems are of particular interest because, by 
tuning the gates, it is possible to change the charge 
on each dot separately. Unfortunately, the aforemen- 
tioned experiments concerned only DQDs coupled to 
nonmagnetic leads. There are several experimental real- 
izations of single quantum dots attached to ferromagnetic 
le^^LS^S^M^SMmiJMS. while experimental 
data on spin-polarized transport through double quan- 
tum dots are lacking. We thus believe that the results 
presented in this paper will be of assistance in discussing 
future experiments. 

The paper is organized as follows: The model of the 
considered system is presented in Sec. II, whereas the 
method employed in calculations is described in Sec. III. 
In Section IV we present and discuss the numerical results 
for double quantum dots coupled in series. In this section 
we also focus on the role of second-order processes in 
transport. In addition, we also analyze transport in the 
Pauli spin blockade regime. In Sec. V we deal with 
DQDs coupled in parallel, whereas in Sec. VI we consider 
transport through T-shaped quantum dots. Finally, the 
conclusions are given in Sec. VII. 



II. MODEL 

The schematic of a double quantum dot coupled to fer- 
romagnetic leads is shown in Fig. [TJ It is assumed that 
the magnetizations of the leads are oriented collinearly, 
so that the system can be either in the parallel or an- 
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FIG. 1: (Color online) Schematic of a double quantum dot 
weakly coupled to external ferromagnetic electrodes. The 
magnetizations of the leads can form either parallel or an- 
tiparallel magnetic configuration. The first (j — 1) and sec- 
ond (j = 2) dots are coupled to each other via the hopping t, 
and to the left (r — L) and right (r = R) leads with the cou- 
pling strength T r j . The double quantum dot is assumed to be 
symmetrically biased. By adjusting the couplings, the system 
can smoothly cross over from the serial to parallel geometry. 

tiparallel magnetic configuration. The Hamiltonian of 
the system is given by 

H = Hl + + #DQD + Ht (1) 

The first two terms describe noninteracting itinerant elec- 
trons in the leads, H r = J2ka e rkacl ka c r ] icr for the left 
(r = L) and right (r = R) lead, where e r ko- is the energy 
of an electron with the wave vector k and spin a in the 
lead r, and c^ k(T (c r kcr) denotes the respective creation 
(annihilation) operator. The double dot is described by 
the Hamiltonian 

j=l,2 a j=l,2 

+ U'J2 ni„n 2a > + t ^2(4 a d 2a + 4 a d la ) ,(2) 

with rij a — d^dju, where (dj a ) is the creation (anni- 
hilation) operator of an electron with spin a in the first 
(j = 1) or second (j = 2) quantum dot, and e 3 - is the 
corresponding single-particle energy. The Coulomb in- 
teraction on the first and second dot is assumed to be 
equal and is described by U, while U' corresponds to the 
inter-dot Coulomb correlation. The last term of Hdqb 
describes the hopping between the two dots with t being 
the hopping parameter. We assume that the hopping pa- 
rameter is large, so that there is a considerable overlap 
of the wave functions of the two dots, leading to the for- 
mation of molecular many-body states, through which 
transport takes place. In addition, we also note that 
an exchange interaction between spins in the two dots 
may lead to the formation of singlet and triplet states*^ 
However, experimentally, this exchange interaction was 
found to be rather small as compared to the other en- 
ergy scales^ therefore, in the following considerations, 
we will neglect it. 

The tunneling processes between the DQD and elec- 
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trodes are described by the Hamiltonian, 

r=L,Rj=l,2 k<T 

where t r j denotes the tunnel matrix elements between 
the rth lead and the jth dot. The coupling of the jth 
dot to the rth lead can be written as T°j = 2n\t r j\ 2 p^, 
where p% is the spin-dependent density of states of lead 
r. By introducing the definition of the spin polarization 
of lead r, p r = (p+ — p~)/(p+ + Pr)i one can write, 

r 5 _) = r rj (i±p r ), with r rj = (r+.+r-.)/2, where r+ 

and r~- describe the coupling of the jth dot to the spin- 
majority and spin-minority electron bands of the lead r, 
respectively. 

Because the double dot Hamiltonian, Eq. is not 
diagonal in the local basis, we perform a unitary trans- 
formation, W HdqdII = Hbqd, to a new basis in which 

#dqd is diagonal, #dqd|x) = £ x \x)- Tne eigenvectors 
\x) are the many-body states of the double quantum 
dot, while the eigenvalues e x denote the corresponding 



energies 
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III. METHOD 

In order to calculate the spin-dependent transport 
properties of a double quantum dot in the sequential 
and cotunneling regimes, we employ the real-time dia- 
grammatic techniquei 43 ' 75 i 76 It consists in a systematic 
perturbation expansion of the reduced density matrix of 
the considered system and the current operator with re- 
spect to the dot-lead coupling strength T. The time evo- 
lution of the reduced density matrix can be visualized 
as a sequence of irreducible self-energy blocks, W xx > , on 
the Keldysh contour, where W xx ' describe transitions be- 
tween the many-body DQD states \\) and \x')- The ele- 
ments W xx > set up a self-energy matrix W. Within the 
matrix notation introduced by Thielmann et al.j^ the 
stationary elements of the reduced density matrix can be 
found from 



(Wp 



si \ 

h 



r<5. 



(4) 



where p st is the vector containing probabilities and the 
matrix W is the modified matrix W so as to include 
the normalization of probabilities. The current flowing 
through the system can be then found from 



2fi 



Tr{WV } 



(5) 



The matrix W 1 is the self-energy matrix with one inter- 
nal vertex resulting from the expansion of the tunneling 
Hamiltonian replaced by the current operator. 

To calculate the transport properties order by order in 
tunneling processes, we expand the self-energy matrices, 
W = W« + W (2 ) + . . . , W 1 = W 1 ^) + W 1 ' 2 ) + . . . , and 



the dot occupations, p st = p st (°) + p st (!) + . . . , respec- 
tively. Then, the first-order (sequential) and the second- 
order (cotunneling) currents are given by 



j(D 



2H 

e 



Tr { W I(l) p st(0) } 



jP) = ^{w^p-M+W^yW}. (6) 

On the other hand, the zeroth and first-order occupation 
probabilities can be found from the following equations 



r<s v 



j x — _L UxXo , 

w<V t(1) +w(V t (°> - 0, (7) 

where W« [w( 2 )j is given by W« [W^ 2 )] with one 

arbitrary row \o replaced by (r, . . . , T) [(0, . . . , 0)] due 

to the normalization of probabilities, Yl x Px = ^™' ' 

As one can see from the above formulas, to deter- 
mine the transport properties it is necessary to calcu- 
late the elements of the corresponding self-energy ma- 
trices. This can be done using the respective diagram- 
matic rules i 43 i 75 ' 76 Although the first-order calculation 
is rather simple, the general analytical formulas for the 
self-energies in the second order are rather complicated 
due to many virtual states through which the cotunneling 
processes can take place. In the Appendix A, as an exam- 

(2) 

pie, we present the contribution coming from I^vu WiV) 
where N is the charge state of the double dot. 

In addition, in the present paper we also analyze the 
zero-frequency current noise of the double quantum dot 
system, S = 2j^ dfc((/(t)/(0) + J(0)/(t)) - 2(/) 2 ). For 
low bias voltages the current noise is dominated by ther- 
mal noise, while for \eV\ > k-&T the noise associated with 
the discrete nature of charge (shot noise) dominates. 80 
Within the real-time diagrammatic technique, a general 
formula for the shot noise has been derived in Ref. [76| 
taking into account the non-Markovian effects. The shot 
noise in the first-order can be found from the correspond- 
ing expression 



- ^Tr 

h 



(\y II(1) + W I(1) P (-1) W I(1) ) p st(0) 



while the cotunneling current noise is given by 



(8) 



5(2) = £_Tr(|w II(2) + W I(2) P ( ~ 1} W I(1) 

+ w I(1) P (_1) W I(2) + W I(1) P (0) W I(1) 
+ W^Q^AWWl p st (°) 

+ [w^ + W'WpHlwiWlp^D^ (9) 

where QW = p st (") ® e T , with e T = (1, . . . , 1). The ob- 
jects P ( ~ 1} and P (0) are given by, W^pt" 1 ) = QW - 1, 
and, WWpM +W( 2 )P(-!) = 1(QW -aW^QW), re- 
spectively, with 1 being the unit vector with row xo set to 
zero. On the other hand, the matrices W 11 ^ 1 ) and W 11 ' 2 ' 
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are the first and second-order self-energy matrices with 
two internal vertices replaced by the current operator, 
while <9W« and dW 1 ^ are partial derivatives of W' ' 
and W ! W with respect to the convergence factor of the 
Laplace transform^ 

By taking into account all the first-order and second- 
order contributions to the self-energies, we are able to 
resolve the transport properties in the full range of the 
bias and gate voltages. The first order of expansion with 
respect to the coupling corresponds to sequential tunnel- 
ing, while the second order is associated with cotunncl- 
ing. Sequential tunneling is allowed if the applied bias is 
larger than the threshold voltage, i.e. when the energy 
provided by the transport voltage is comparable with the 
charging energy. Otherwise, the system is in the Coulomb 
blockade regime where sequential tunneling is exponen- 
tially suppressed and the current flows due to cotunncl- 
ing, which involves correlated tunneling through virtual 
states of the system i 77 ' 78 ' 79 

Among different cotunncling processes one can gen- 
erally distinguish two types of processes: non-spin-flip 
ones, which do not affect the DQD state, and the spin- 
flip ones. Furthermore, one can also have double-barrier 
cotunncling, which contributes directly to the current, 
and single-barrier cotunneling, which affects the double 
dot occupations and, thus, indirectly, the current. 



IV. DOUBLE QUANTUM DOTS COUPLED IN 
SERIES 

In this section we present and discuss numerical re- 
sults on double quantum dots coupled in scries. In order 
to realize this geometry we set IY2 = Tri = and as- 
sume Tli = Tri2 = r/2, see Fig. [T] Furthermore, in 
the following we will also distinguish between symmetric 
and asymmetric DQDs. In the former case, the level po- 
sition of each dot is the same, E\ = £2, while in the latter 
case the levels are detuned, £\ ^ £2- In addition, in this 
section we will emphasize the role of second-order pro- 
cesses in transport. In order to ascribe observed features 
to respective tunneling processes, we will therefore also 
present results obtained within the sequential tunneling 
approximation, i.e. when considering only the first-order 
tunneling processes. 



A. Symmetric double quantum dots 

The differential conductance as a function of the bias 
voltage and the position of the dot levels for the parallel 
and antiparallel configurations in shown in Fig. Ufa) and 
(b). Because experimentally the position of the dots' lev- 
els can be changed upon applying a gate voltage, Fig. [3] 
effectively shows the bias and gate voltage dependence of 
the conductance. First of all, it can be seen that the dif- 
ferential conductance displays a characteristic Coulomb 
diamond structure. The diamonds in the low bias volt- 



1 




-1 




-4 



-5 

-8 -4 4 8 




-8-4 4 




-8 -4 4 8 

V (mV) 



FIG. 2: (Color online) The differential conductance in the par- 
allel Gp (a) and antiparallel Gap (b) magnetic configurations 
as a function of the bias voltage V and the position of the dots' 
levels e = ei — E2 for double quantum dots coupled in series. 
The parameters are: fcp,T = 0.15 meV, U — 2 meV, U' = 1 
meV, t = 0.25 meV, r L2 = T R i = 0, lYi = r R2 = T/2, with 
r = 0.1 meV, and p = 0.5. Because the double quantum dot 
is symmetric, e\ = £2, there is no asymmetry associated with 
the bias reversal. For comparison in part (c) we also show 
the density plot of the differential conductance in the parallel 
configuration calculated using only the first-order tunneling 
processes, Gp . 



age regime correspond to the Coulomb blockade regime 
where the sequential current is exponentially suppressed 
and the current flows due to cotunneling. With increas- 
ing the bias voltage the excited states start participating 
in transport, which leads to additional lines in the dif- 
ferential conductance. The features described above are 
rather associated with the energy spectrum and charge 
states of the DQD than with the ferromagnetism of the 
leads, therefore they are present in both magnetic con- 
figurations. The spin dependence stemming from ferro- 
magnetic electrodes gives rise to a difference of the con- 
ductance in the parallel and antiparallel configurations, 
see Fig. [^a) and (b). The conductance in the parallel 
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FIG. 3: (Color online) The total (first and second order) TMR 
(a) and the first-order TMR (b) as a function of the bias 
voltage V and the level position e for parameters the same as 
in Fig. [2] The two figures are plotted in the same scale. 



configuration is generally larger than that in the antipar- 
allel configuration. This is due to the spin asymmetry 
in the couplings when the leads' magnetizations are an- 
tiparallel. As in the parallel configuration the majority 
(minority) electrons of, let us say, left lead tunnel to the 
majority (minority) electron band of the right lead, in 
the antiparallel configuration the situation is reversed - 
the majority (minority) electrons of the left lead tun- 
nel to the minority (majority) electron band of the right 
lead. This leads to the difference of the conductance in 
both magnetic configurations, see Fig. [2j and to the cor- 
responding tunnel magnetoresistance effect. In addition, 
in Fig. [^c) we also present the bias and gate voltage de- 
pendence of the differential conductance for the parallel 
configuration calculated within the sequential tunneling 
approximation, Gp\ By comparing Fig. dja) and (c), 
one can see that cotunneling leads to broadening of the 
resonance peaks, lowering the magnitude of the differen- 
tial conductance, and gives rise to finite conductance in 
the Coulomb blockade regimes. 

The bias and gate voltage dependence of the TMR is 
shown in Fig.[3ja), while Fig.[3jb) presents the first-order 
TMR, TMR (1) . Intuitively, the TMR is a measure of the 
system's transport properties change when the magnetic 
configuration switches from the parallel to the antiparal- 
lel one. It is defined as2L2&^2 



TMR = 



/At 



(10) 




where Ip (Zap) is the current flowing through the sys- 



FIG. 4: (Color online) The linear conductance in the parallel 
(solid line) and antiparallel (dashed line) magnetic configura- 
tions (a) and the linear TMR (b) as a function of the level 
position e for parameters the same as in Fig. [2] For com- 
parison, the dotted curves show the linear conductance and 
the TMR calculated taking into account only the first-order 
processes. 



tern in the parallel (antiparallel) configuration. By com- 
parison with Fig. [21 one can easily identify the different 
transport regimes. Furthermore, by comparing the total 
and sequential TMR, one can immediately see that co- 
tunneling modifies the TMR mainly in the low bias volt- 
age regime, i.e. in the Coulomb blockade regime and in 
the regime where sequential tunneling is suppressed due 
to the absence of levels in the energy window provided 
by transport voltage. In order to discuss and see more 
clearly the behavior of the TMR on applied voltages, in 
Fig. 0] we show the linear response conductance and re- 
sulting TMR, while in Fig. [5] we display the transport 
properties in the nonlinear response regime. 

The linear conductance in both magnetic configura- 
tions is shown in Fig. 0Ja). When lowering the position 
of the DQD levels the conductance displays four reso- 
nance peaks associated with subsequent occupation of 
the corresponding charge states. The dotted curves in 
Fig. 2] present results obtained in the first-order approx- 
imation. It can be seen that off resonance the current is 
mainly dominated by cotunneling. On the other hand, 
on resonance the sequential processes dominate current, 
while cotunneling only slightly affects the conductance, 
leading to renormalization of the DQD levels and, thus, 
slightly shifting the position of the conductance peaks, 
see Fig. Ufa). Interestingly, the second-order processes 
have a rather large impact on the linear TMR shown in 
Fig. IUb). First of all, the linear response TMR exhibits 
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a strong dependence on the gate voltage, i.e. on the 
number of electrons in the double quantum dot. When 
the DQD is either empty or fully occupied, the TMR 
is given by the Julliere value, TMR = 2p 2 /(I — p 2 )£L 
This is due to the fact that in those transport regimes 
the current is driven by elastic cotunneling processes 
which do not affect the DQD state in any way. Thus, 
as far as the TMR is concerned, the system behaves as a 
single ferromagnetic tunnel junction^ However, in the 
Coulomb blockade regions, the TMR becomes generally 
suppressed due to the presence of inelastic cotunneling 
which introduces spin-flip processes in the system. The 
dotted line in Fig.fJJb) shows the TMR calculated within 
the sequential tunneling approximation - it is given by 
TMR = p 2 /(l - p 2 )M& At this point, we would like 
to note that the results for the linear conductance cal- 
culated by including only the first-order processes are 
rather reliable in the whole range of e. On the other 
hand, the results for the linear response TMR are com- 
parable to those obtained within the sequential tunneling 
approximation only on resonance, where sequential pro- 
cesses dominate, while off resonance they are completely 
unreliable, see Fig. [¥] Therefore, in order to properly 
analyze the dependence of the TMR in the full range 
of parameters one has to take into account cotunneling 
processes. 

The bias voltage dependence of the current, differen- 
tial conductance, tunnel magnetoresistance and the Fano 
factor is shown in Fig. [5] for e = —2 meV and e = 2 
meV. The first situation, e — — 2 meV, corresponds to 
the case when the ground state of the double quantum 
dot is doubly occupied, while the second one, e — 2 meV, 
corresponds to the case of empty double dot. In all cases, 
due to the spin asymmetry of tunneling processes, the 
current in the parallel configuration is larger than the 
current in the antiparallel configuration. In addition, the 
i-v curves display characteristic Coulomb steps. With 
increasing the bias voltage more and more charge states 
start to participate in transport which gives rise to the 
corresponding steps in the current-voltage characteris- 
tics, and peaks in the differential conductance. Further- 
more, the single-electron charging effects also lead to the 
oscillatory-like behavior of the TMR effect, see Fig. [5jc) 
and (g). In the Coulomb blockade regime, e = — 2 meV, 
the TMR is much suppressed as compared to the case of 
e = 2 meV. This is due to the presence of spin-flip cotun- 
neling processes in the case of the doubly occupied DQD, 
as already discussed in association with Fig. |U However, 
when increasing the bias voltage the TMR increases and, 
at the threshold for sequential tunneling, reaches a local 
maximum, see Fig.[5][c). This effect is due to the nonequi- 
librium spin accumulation which is induced in the DQD 
system with increasing the transport voltage. (Similar 
effect can be observed when the double dot is singly oc- 
cupied, i.e. when e — —0.5 meV.) On the other hand, 
when the DQD is empty, the TMR at low bias voltage 
acquires the Julliere valued and then, with increasing 
V, becomes generally suppressed, see Fig. [SJg). 




FIG. 5: (Color online) The current in units of Io ~ eT/h (a,e), 
differential conductance (b,f) and the Fano factor (d,h) in 
the parallel (solid line) and antiparallel configuration (dashed 
line) as well as the TMR (c,g) as a function of the bias voltage 
for different values of the level position e as indicated in the 
figure. £ = — 2 meV (a-d) corresponds to the case when the 
ground state of the DQD is doubly occupied, while for e = 2 
meV (e-h) the double dot is empty. The parameters are the 
same as in Fig. [2] The dotted curves show the results obtained 
in the sequential tunneling approximation. 



In addition, in Fig. O we also show the Fano factor, 
F = Sy(2e|I|), calculated in both magnetic configura- 
tions and for two values of e, as indicated in the fig- 
ure. First of all, one can see that the Fano factor be- 
comes divergent at low bias voltages. This is due to 
the thermal noise which dominates the current noise as 
V — > 0, while 1^0, leading to F — > cc<££ Further- 
more, it can be seen that the Fano factor in the Coulomb 
blockade regime is slightly larger in the parallel config- 
uration than in the antiparallel one. This is associated 
with an additional positive cross-correlations, which con- 
tribute to the current noise, due to the ferromagnetism 
of the electrodes. In the Coulomb blockade regime, there 
is a larger asymmetry between the cotunneling processes 
through the majority and minority spin channels in the 
parallel configuration than in the antiparallel one. This 
effectively increases the current fluctuations in the par- 
allel configuration and leads to the corresponding dif- 
ference in the Fano factors, see Fig. E^d). On the other 
hand, in the sequential tunneling regime, the information 
about the spin asymmetry of tunneling processes between 
the DQD and the corresponding lead is rather contained 
in the nonequilibrium spin accumulation induced in the 
double dot. This spin accumulation is larger in the an- 
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V(mV) 

FIG. 6: (Color online) The total Fano factor in the parallel 
(a) and antiparallel (b) magnetic configurations as a function 
of the bias voltage V and the position of the dots' levels e 
for parameters the same as in Fig. [2] Part (c) presents the 
Fano factor in the parallel configuration calculated in the first- 
order approximation. Because the Fano factor is divergent 
for \eV\ < kuT, this transport regime is denoted with a thick 
black line. 



tiparallel configuration, yielding generally an enhanced 
Fano factor in the antiparallel configuration as compared 
to the parallel one. Moreover, we also note that in the 
Coulomb blockade regime, e = —2 meV, the shot noise 
becomes super-Poissonian (F > 1) and drops to the sub- 
Poissonian value (F < 1) at the threshold for sequential 
tunneling, see Fig.[5jd). This super-Poissonian shot noise 
is associated with bunching of inelastic cotunncling pro- 
cesses through the system.! 54 ' 81 On the other hand, in 
the case of empty DQD, e = 2 mcV, in the cotunnel- 
ing regime the current flows due to elastic second-order 
processes which obey the Poissonian statistics. In this 
case the shot noise is given by S = 2e\I\ and the Fano 
factor is simply equal to unity, see Fig.[5Jh), irrespective 



of magnetic configuration of the system. As concerns 
the sequential tunneling regime, the shot noise is then 
generally sub-Poissonian, irrespective of the position of 
the DQD level e. This indicates the role of correlations 
in electronic transport, in particular, the Coulomb cor- 
relations and charge conservation^ Our results are in 
qualitative agreement with experimental data on tunnel- 
ing through vertically coupled self-assembled InAs quan- 
tum dots where also super-Poissonian shot noise has been 
founds 

Finally, to make the present analysis self-contained, we 
also show the bias and gate voltage dependence of the to- 
tal Fano factor for the parallel and antiparallel magnetic 
configuration, see Fig.^a) and (b), respectively. In addi- 
tion, for comparison, the Fano factor calculated by taking 
into account only first-order tunneling processes is de- 
picted in Fig. [6](c). The black lines around the zero bias, 
\eV\ < ksT, mark the transport regime where the Fano 
factor is divergent due to finite thermal noise. The differ- 
ent behavior of the shot noise is now clearly visible. In the 
Coulomb blockade regime the noise is super-Poissonian, 
in the cotunneling regime when the DQD is either empty 
of fully occupied (current is mediated by elastic cotun- 
neling) the shot noise is Poissonian, and in the sequential 
tunneling regime the noise drops to sub-Poissonian value. 



B. Asymmetric double quantum dots: Pauli spin 
blockade 

By applying a gate voltage to each quantum dot, it 
is possible to tune the dot levels separately. So far, we 
have considered the case of symmetric DQD, i.e. when 
£\ = £2- In this situation the transport characteristics 
were symmetric with respect to the bias reversal. This 
is however not the case for s\ =/= £2, where the current 
becomes asymmetric with respect to the bias reversal, 
leading to the Pauli spin blockade and negative differen- 
tial conductance, as observed experimentally 6,8 In the 
following, assuming realistic parameters of the double 
quantum dot system, we analyze transport properties in 
the regime where the Pauli spin blockade effects are vis- 
ible. 

The mechanism leading to the spin blockade was theo- 
retically discussed by Fransson et al^ and Muralidharan 
et al^ These considerations were however restricted only 
to the first-order tunneling processes, which dominate the 
current out of the blockade regime. As shown experimen- 
tally by Ono et al.^ in the spin blockade regime there is a 
finite leakage current, which cannot be explained within 
the sequential tunneling approximation. To explain the 
existence of this leakage, it was proposed that non-zero 
current in the Pauli blockade is associated with spin-flip 
processes induced by hyperfine interaction^ In the fol- 
lowing, we show that the leakage current results just from 
the interplay between different intrinsic tunneling pro- 
cesses driving the current. The hyperfine interaction, or 
coupling to a phonon bath, may, of course, increase the 
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FIG. 7: (Color online) The current (a) and Fano factor (c) in 
the parallel (solid) and antiparallel (dashed) magnetic config- 
urations as well as the TMR (b) as a function of the bias 
voltage for DQD coupled in series. The parameters are: 
k B T = 0.05 meV, ei = -1 meV, e 2 = -2 meV, A = 3 meV, 
U = 2 meV, U' = 1 meV, t = 0.05 meV, where A is the level 
spacing in each dot. The other parameters are the same as in 
Fig. [2] The inset in part (a) displays the current in the paral- 
lel magnetic configuration in the Pauli spin blockade regime. 
The dotted lines correspond to the first-order calculation. 



leakage, but is not necessary for the observation of a finite 
current in the Pauli spin blockade regime. 

When the DQDs levels are detuned, £2 < £1, and once 
the first and second dot become singly occupied, the cur- 
rent may be suppressed in some range of the bias voltage. 
This is associated with the full occupation of two-electron 
triplet states of the DQD^ The Pauli spin blockade may 
be lifted when the applied bias voltage admits another 
(exited) charge states to participate in transport. On the 
other hand, when the voltage is reversed, it is energeti- 
cally allowed that the electron from the first dot tunnels 
to the left lead and then another electron from the right 
lead enters the DQD, lifting the Pauli spin blockade. 



This can be seen in Fig. [7] where we present the bias 
voltage dependence of the current, Fano factor and the 
TMR effect. In order to make the calculations more re- 
alistic and to allow for excited states of the system, we 
have now taken into account four different orbital levels, 
two in each dot. Such system is described by straight- 
forward extension of the hamiltonian, Eq. where the 
level spacing in the dots is described by the parameter A. 
First of all, one can see that the transport characteristics 
are asymmetric with respect to the bias reversal. This is 
associated with the fact that by detuning the DQDs lev- 
els the symmetry of transition rates between the left and 
right leads has been broken. Furthermore, for positive 
bias voltage the current is suppressed in a broad range of 
the bias voltage (0 < eV < U) due to the full occupation 
of the triplet states, while for negative bias the Pauli spin 
blockade is lifted. In addition, the blockade can also be 
lifted when the positive bias voltage is increased further, 
eV > U, so that tunneling through the second level of 
the second dot is allowed, see Fig. [TJa). 

In the inset of Fig. EJa) we show the current in the 
parallel configuration just in the spin blockade regime. 
For comparison, we also display the current calculated 
within the sequential tunneling approximation. One can 
see that the sequential current is suppressed as compared 
to the total current (calculated taking into account co- 
tunneling processes). The suppression of the sequential 
current is generally governed by the ratio 2t/|e 2 — ei\M> 
When 2t/\s2— £i| <C 1, the occupation of the two-electron 
triplet states approaches unity and the first-order current 
becomes fully blocked. However, the second-order pro- 
cesses are still allowed, leading to a finite leakage current 
in the spin blockade regime. One can distinguish different 
contributions coming from cotunneling: 

(i) The double-barrier elastic second-order processes 
which contribute directly to the current. 

(ii) The double-barrier spin-flip cotunneling which con- 
tributes to the current and reduces the occupation of the 
triplet state, this way opening the DQD for the sequential 
tunneling. 

(iii) The single-barrier spin-flip cotunneling which does 
not contribute directly to the current but flips the spin 
in the DQD and allows the sequential processes to occur. 

As a consequence, in the Pauli spin blockade regime 
the current flows due to cotunneling and spin-flip 
cotunncling-assisted sequential tunneling. The key role 
is played by the spin-flip processes which decrease the 
occupation of the two-electron triplet states and lead to 
a finite occupation of singlet states. The importance of 
the spin-flip processes in explanation of the finite leakage 
current has been invoked in Ref. [171 ] . The leakage was 
there associated with spin-flip processes induced by hy- 
pcrfine interaction through the Overhauser effect. Here, 
we show that the leakage results just from the interplay 
between different tunneling processes, i.e. just from the 
pure nature of tunneling processes. We also note that the 
effects of cotunneling are rather independent of the mate- 
rial from which the dots are built, unlike the Overhauser 
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FIG. 8: The total (sequential plus cotunneling) current (solid 
line) and the sequential current (dotted line) in the paral- 
lel configuration as a function of the hopping between the 
two dots calculated in the middle of the Pauli spin blockade 
regime, V — 1 mV. The parameters are the same as in Fig. [7] 
The dependence of the current on 2t/\e2 — ei| in the antipar- 
allel configuration is qualitatively similar. 



fields Thus, for example, in GaAs DQDs the leakage 
current in the Pauli spin blockade would be associated 
with both cotunneling and hyperfine contributions, while 
in carbon nanotube DQDs it would be mainly due to co- 
tunneling. Furthermore, as shown in experiments by Ono 
and Tarucha, 83 who studied the dependence of the leak- 
age current on the applied magnetic field, the hyperfine 
interaction starts to play a role at certain finite in-plane 
magnetic fields, at which the nuclei become polarized, 
leading to a sudden jump of the leakage current as a 
function of magnetic field. This suggests that the main 
contribution to the current in the spin blockade regime 
in the absence of magnetic field may come from cotun- 
neling. In the weak coupling regime, typical values of the 
dot-lead coupling strength T are of the order of /ieV, 84 
which gives the leakage current of the order of 10 _3 Jo « 1 
pA (see Fig. [8]), in agreement with experimental results^ 3 - 
Finally, we notice that to explain the jump in the leakage 
current when sweeping magnetic field, besides cotunnel- 
ing, one also has to include the hyperfine interaction. 
Nevertheless, this goes beyond the scope of the present 
paper. 

In Fig. [8] we present the logarithmic dependence of the 
sequential and total currents in the parallel configura- 
tion on the hopping parameter t. It is clearly visible that 
for 2i/|e 2 — £i | <C 1, the sequential current is smaller 
by two orders of magnitude than the total current. On 
the other hand, when 2t/\e2 — £i\ ~ 1, the sequential 
current becomes of the same order as the total current, 
but still, noticeably, the total current is larger than that 
calculated using only the first-order tunneling processes. 
The saturation of the two currents for 2i/|e 2 — £i| -C 1, 
see Fig. HI is associated with a finite temperature. How- 
ever, because sequential tunneling depends exponentially 
on temperature while cotunneling only algebraically, at 
lower T the difference between the two currents would be 



even more pronounced. The dependence of the current 
in the antiparallel configuration is qualitatively similar to 
the one shown in Fig. [8l We also note that generally the 
above mentioned mechanism responsible for the leakage 
current does not depend on whether the leads are made 
of ferromagnetic or nonmagnetic material. 

The TMR as a function of the bias voltage is displayed 
in Fig. [Tfb). One can see that the TMR calculated us- 
ing the first-order processes and that calculated taking 
into account cotunneling are generally similar, except for 
the Pauli spin blockade regime. In this transport regime 
the sequential TMR is much overestimated. To under- 
stand this behavior we note again that the blockade re- 
sults from the full occupation of DQD's triplet states. 
The probability is equally distributed between the three 
components of the triplet in the parallel configuration. 
However, in the antiparallel configuration, it turns out 
that the charge is mainly accumulated in the S z = 1 
component of the triplet. As a consequence, in the par- 
allel configuration all the three components of the triplet 
participate in transport, while in the antiparallel config- 
uration only one. Because the sequential current in the 
blockade regime is mainly associated with thermal fluc- 
tuations, the difference in the number of states relevant 
for transport leads to large sequential TMR in blockade 
regime. This difference also gives rise to the correspond- 
ing difference between Fano factors calculated in the se- 
quential tunneling approximation, see the dotted curves 
in Fig.[7Jc). The shot noise is then super-Poissonian and 
larger in the parallel configuration. However, as pointed 
above, the leakage current in the Pauli spin blockade 
regime is due to cotunneling and cotunneling-assisted 
sequential tunneling processes. The total current flows 
not only due to thermal fluctuations but due to corre- 
lated tunneling through virtual states of the DQD. This 
fact generally decreases the current fluctuations and the 
difference between the two magnetic configurations. As 
a result, the total TMR and Fano factor become sup- 
pressed as compared to the sequential tunneling results, 
see Fig. [7Jb) and (c), respectively. In addition, we also 
notice that, irrespective of the magnetic configuration of 
the system, the shot noise is super-Poissonian in the spin 
blockade regime and drops to sub-Poissonian value out of 
the blockade regime, see Fig.JT^c). The super-Poissonian 
shot noise is generally due to bunching of spin-flip co- 
tunneling and cotunneling-assisted sequential tunneling 
processes. 



V. DOUBLE QUANTUM DOTS COUPLED IN 
PARALLEL 

In this section we analyze the spin-polarized transport 
through double quantum dots coupled in parallel. This 
geometry can be realized by setting T T j = T/2, for r — 
L, R and j = 1, 2, see Fig. [T] We note that in the case of 
negligible hopping between the two quantum dots, t — > 
0, the behavior of the DQD system resembles that of a 
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single multi-level quantum dot. The problem of spin- 
dependent transport through multi-level quantum dots 
has been addressed very recently in Ref . [5J] and will not 
be considered here. 

In Fig. [9] we present the density plots of the differential 
conductance and Fano factor in the parallel configuration 
as well as of the TMR effect as a function of the bias 
and gate voltages. The behavior of the conductance and 
the shot noise in the antiparallel configuration is qualita- 
tively similar as in the parallel configuration, therefore, it 
is not shown here. The information about the difference 
in transport in the two magnetic alignments is contained 
in the TMR, whose magnitude reflects the asymmetry 
in tunneling processes when the leads are parallel or an- 
tiparallel to each other. 

The differential conductance, shown in Fig. [Ha), dis- 
plays characteristic Coulomb diamonds. Because the in- 
ternal energy structure of the DQD is generally the same 
as in the case of DQD coupled in series, the differential 
conductances are qualitatively similar, see Figs. [2a) and 
UK a). The main difference is that in the case of parallel 
DQDs the level of each dot is coupled both to the left 
and right leads, which generally results in an enhanced 
conductance for parallel DQDs as compared to DQDs 
coupled in series. Similarly the TMR shown in Fig.[9jb), 
its bias and gate voltage dependence is only slightly mod- 
ified as compared to Fig. [3ja). Again, the linear TMR 
exhibits a strong dependence on the occupation number 
of the DQD, with TMR given by the Julliere formula 
for empty and fully occupied double quantum dot, and 
is much suppressed in other transport regimes due to 
spin-flip cotunneling. In the Coulomb blockade regime, 
the transport properties are mainly conditioned by the 
spin-flip cotunneling, which gives a dominant contribu- 
tion to the current. Furthermore, the bunching of spin- 
flip cotunneling in the Coulomb blockade regime leads to 
super-Poissonian shot noise, see Fig. [5Jc). The super- 
Poissonian shot noise is however more pronounced in 
Coulomb blockade regimes with an odd number of elec- 
trons in the double dot. On the other hand, for voltages 
above the threshold for sequential tunneling, the shot 
noise becomes sub-Poissonian, with the Fano factor ap- 
proaching 1/2. 

Summing up, we notice that generally in the weak 
coupling regime the effect of different geometries of the 
double dot system leads only to qualitative difference in 
transport properties. The main difference is in the mag- 
nitude of conductance - for parallel DQDs it is roughly 
two times larger than for serial DQDs. This is contrary 
to the strong coupling regime where, for example, in the 
parallel geometry the orbital Kondo phenomenon arises 
due to the interference effects, while for DQDs coupled in 
series the orbital Kondo effect is destroyed^ Finally, we 
also note that in the case of detuned levels, |e2 — £i | 3> 0, 
transport properties of double quantum dots coupled in 
parallel are qualitatively similar to those of single multi- 
level quantum dots^ 
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FIG. 9: (Color online) The differential conductance (a) and 
Fano factor (c) in the parallel magnetic configuration, and 
the TMR (b) as a function of the bias voltage V and the 
position of the dots' levels e = ej = £2 for double quantum 
dots coupled in parallel. The parameters are the same as in 
Fig. [2] with r rj = T/2, for r = L,R, j = 1,2, and T = 0.1 
meV. The transport regime where the Fano factor is divergent 
is marked with a thick black line. 



VI. T-SHAPED DOUBLE QUANTUM DOTS 

An interesting situation occurs when only one of the 
two quantum dots is coupled to external leads, while the 
other one is decoupled. In such T-shaped systems trans- 
port takes place through the molecular states of the dou- 
ble dot, although only the first dot is directly coupled 
to the leads. This may lead to new transport behavior, 
especially to a large super-Poissonian shot noise and an 
enhanced TMR, as we show in the sequel. 

In Fig. [10] we present the current, differential conduc- 
tance, the TMR and the Fano factor as a function of the 
bias voltage for the case when the first dot is coupled 
to the leads while the second dot is decoupled. Spin- 
dependent transport through single quantum dots have 
already been extensively studied. However, the hopping 
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FIG. 10: (Color online) The bias dependence of the current 
and differential conductance (a), the TMR (b), and the Fano 
factor (c) in the parallel (solid line) and antiparallel (dashed 
line) configurations for T-shaped double quantum dots. The 
parameters are: r r i = T/2, and F r 2 = for r = L,R, with 
r = 0.1 meV, and e x = 1 meV, e 2 = meV, k B T = 0.18 
meV, U = 2 meV, U' = 1 meV, t = 0.2 meV, and p = 0.5. 
The dotted curves present the TMR (b) and the Fano factor 
in the parallel configuration (c) calculated for F r 2 = r r i. 



t between the two dots may modify the transport prop- 
erties significantly, giving rise to novel behavior. First 
of all, we note that due to the finite t, transport takes 
place through molecular states of the double dot. This 
in turn leads to additional steps in the i-v curves and 
peaks in the differential conductance, see Fig. fTUT a). as 
compared to tunneling through one single-level quantum 
dots, where for spin-degenerate level the current exhibits 
only two steps. The difference between the currents in the 
parallel and antiparallel configurations gives rise to the 
TMR which is shown in Fig.QUb). The TMR displays an 
oscillatory-like behavior with increasing the bias voltage. 
In addition, at the threshold for sequential tunneling, the 
TMR becomes increased above the Julliere TMR^i This 
result is rather counterintuitive, as the Julliere TMR is 



FIG. 11: (Color online) The bias dependence of the current, 
differential conductance (a), TMR (b) and the Fano factor 
(c) in the parallel (solid line) and antiparallel (dashed line) 
configurations. The parameters are the same as in Fig. [10] 
with £2 = —4 meV. The dotted curves present the TMR (b) 
and the Fano factor in the parallel configuration (c) calculated 

for r r2 = r ra . 



characteristic of a single ferromagnetic tunnel junction 
and one would expect that an additional object (quan- 
tum dot) between the leads should decrease the TMR. 
This was shown for single quantum dots with spin degen- 
erate levels and symmetric couplings to the leads where 
the TMR was found to take at most the Julliere valued 
In principle, the TMR can be enhanced above the Jul- 
liere value when only one spin component of the dot takes 
part in transport. Such spin selection may be achieved 
upon applying an external magnetic field or due to finite 
exchange interaction between spins in the doti 46 i 86 In the 
case of T-shaped DQDs, we find that the enhancement of 
the TMR is associated with increased occupation of the 
second dot in the antiparallel configuration, as compared 
to the parallel. Because the second dot is decoupled from 
the leads, this increases the difference between the cur- 
rent in the two configurations, yielding the TMR larger 
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than the Julliere value. In addition, finite occupation of 
the decoupled dot also leads to large current fluctuations, 
as the rate for tunneling between the second dot is much 
slower than that for tunneling between the first dot and 
the leads. This in turn gives rise to super-Poissonian shot 
noise, as displayed in Fig. riffl e). The super-Poissonian 
shot noise is present in both magnetic configurations. It 
is also slightly larger in the parallel configuration than in 
the antiparallel one, which is due to an additional contri- 
bution to the noise coming from spin-dependent tunnel- 
ing between the dots and the leads^ We note that the 
super-Poissonian shot noise in T-shaped dots has also 
been reported in the case of nonmagnetic leads^i The 
two dots were however rather weakly coupled to each 
other so that no molecular states were formed, contrary 
to the case considered here. 

The enhancement of the TMR above the Julliere value 
and the large super-Poissonian shot noise are present 
approximately at the threshold for sequential tunneling. 
In this transport regime the current is mainly mediated 
through the charge states of the decoupled dot. The 
aforementioned effects should therefore become washed 
out if there were a finite coupling between the second dot 
and the leads. This is shown in Fig. [Tfflb) and (c) where 
we depict the TMR and Fano factor in the parallel con- 
figuration calculated for r r 2 = T r i (see the dotted curves 
in Fig. [TDD. The TMR larger than the Julliere TMR 
accompanied by super-Poissonian shot noise is thus an 
indication that transport takes place through quantum 
dot coupled to another dot, which is not connected to 
the leads directly. 

Figure [Til displays the bias dependence of the current, 
TMR and the Fano factor calculated in the case when 
at equilibrium the second dot is fully occupied while the 
first dot is empty. At low bias voltage, the system is in 
the Coulomb blockade and the double dot is occupied 
with two electrons in the second dot with unit probabil- 
ity. The current flows then only due to elastic cotunneling 
processes. This gives rise to the Julliere TMR and Pois- 
sonian shot noise, see Fig. [TTT b) and (c). With increasing 
the bias voltage, close to the threshold for sequential tun- 
neling, the occupation of the second dot is slightly low- 
ered at the cost of a finite occupation of another states. 
This opens the system for the cotunneling-assisted se- 
quential and sequential tunneling processes. Neverthe- 
less, because the double dot is still mainly occupied by 
two electrons on the decoupled dot, this gives rise to ex- 
tremely large current fluctuations and Fano factors, see 
Fig. [TlTc). This super-Poissonian shot noise can be con- 
siderably reduced once the second dot becomes coupled 
to the leads, see the dotted curve in Fig. [TTT ch 

Finally, in Fig. rT2] we present the differential conduc- 
tance in the parallel configuration and the TMR as a 
function of the bias voltage and the position of the sec- 
ond dot level. The Coulomb blockade regimes are clearly 
visible in Fig. [T27 a). whereas the different behavior of the 
TMR depending on the transport regime is presented in 

Fig. Hub). 
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FIG. 12: (Color online) The differential conductance (a) in 
the parallel magnetic configuration, and the TMR (b) as a 
function of the bias voltage V and the position of the second 
dot level e = £2 for T-shaped double quantum dots. The 
parameters are the same as in Fig. 1101 



VII. CONCLUSIONS 

We have analyzed the spin-dependent transport 
through systems built of two strongly coupled quantum 
dots which are weakly connected to external ferromag- 
netic leads. The considerations were based on the real- 
time diagrammatic technique which allowed us to deter- 
mine transport properties in the sequential and cotunnel- 
ing regimes in a fully systematic way. In particular, we 
have analyzed the current, differential conductance, shot 
noise and the TMR for different geometries of the double 
dot system, including the serial and parallel couplings, as 
well as the T-shaped systems. We have also discussed the 
main differences in transport characteristics correspond- 
ing to different DQDs geometries, which may be helpful 
in determining the geometry of coupled quantum dots in 
experiments. 

In the case of double quantum dots coupled in series, 
we have found an interesting dependence of the TMR 
on the occupation number of the DQD. Furthermore, 
the super-Poissonian shot noise in the Coulomb blockade 
regimes have been observed. On the other hand, when 
the levels of the DQD were detuned, transport character- 
istics revealed the Pauli spin blockade effects. We have 
shown that the leakage current, observed experimentally 
in the blockade regime, 6 results from the interplay of co- 
tunneling processes which flip the spin in the DQD and 
make the sequential tunneling possible. Thus, the current 
in the spin blockade flows due to cotunneling and spin- 
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flip cotunncling-assisted sequential tunneling processes. 
This mechanism is associated with pure nature of tunnel- 
ing processes and is thus relevant for both nonmagnetic 
and ferromagnetic leads. In addition, we have also shown 
that the shot noise in the Pauli spin blockade regime be- 
comes super-Poissonian, while outside the blockade it is 
sub-Poissonian. 

For double quantum dots coupled in parallel, transport 
characteristics were found to be qualitatively similar to 
those in the case of DQDs coupled in series. The main 
difference is in the magnitude of conductance - for par- 
allel DQDs it is roughly two times larger than for serial 
DQDs. 

In addition, we have also analyzed the case when the 
first quantum dot is coupled to the leads while the sec- 
ond one is completely decoupled. In such T-shaped sys- 
tems we have found a large super-Poissonian shot noise at 
the threshold for sequential tunneling and the TMR en- 
hanced above the Julliere value. These effects are associ- 
ated with an increased occupation of the decoupled quan- 
tum dot, and become washed out once there is a finite 
coupling between the second dot and the leads. Thus, the 
enhanced TMR together with large super-Poissonian shot 
noise may be an indication that transport takes place 
through quantum dot which is side-coupled to another 



dot disconnected from the leads. 
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APPENDIX A: DETAILS OF NUMERICAL 
CALCULATIONS 

In order to calculate the transport properties in the 
sequential and cotunncling regimes, it is necessary to 
find the elements of the respective first-order and second- 
order self-energy matrices. This can be done using the 
respective diagrammatic rules i 43 ' 75 ! 76 Here, as an exam- 
pie, we present the contribution coming from W^/ N \ x '/m 
where N is the charge state of the double dot. It is given 

by 
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where 7 ^ ± (e) = r^-/ ± (e - Hr)/{2ir), with /+ being the 
Fermi distribution function, /~ = 1 — / + , fi r denoting 
the electrochemical potential of the lead r, and 

Bgffi (e, e') = J durtf ( U - e<) . 

We note that having found all the first-order and second- 
order self-energy matrices, we are in principle able to 



calculate transport through arbitrary number of differ- 
ent orbital levels coupled to each other and to external 
leads in an arbitrary way. As far as numerical details 
are concerned, for systems consisting of larger number of 
orbital levels, in calculations we make use of the block 
structure of the initial Hamiltonian in the charge space, 
and perform unitary transformation of the Hamiltonian 
and the relevant local operators for each block separately. 
In addition, we also store the respective matrix elements 
in blocks labelled by charge quantum numbers. 



14 



* Electronic address: weymann@amu.edu.pl 

1 F. R. Waugh, M. J. Berry, D. J. Mar, R. M. Westervelt, 
K. L. Campman, and A. C. Gossard, Phys. Rev. Lett. 75, 
705 (1995). 

2 J. J. Palacios and P. Hawrylak, Phys. Rev. B 51, 1769 
(1995). 

3 R. H. Blick, R. J. Haug, J. Weis, D. Pfannkuche, K. v. 
Klitzing, and K. Eberl , Phys. Rev. B 53, 7899 (1996). 

4 H. Imamura, H. Aoki, and P. A. Maksym, Phys. Rev. B 
57, R4257 (1998). 

5 R. Ziegler, C. Bruder, and H. Schoeller, Phys. Rev. B 62, 
1961 (2000). 

6 K. Ono, D. G. Austing, Y. Tokura, and S. Tarucha, Science 
297, 1313 (2002). 

7 W. G. van der Wiel, S. De Franceschi, J. M. Elzerman, T. 
Fujisawa, S. Tarucha, and L. P. Kouwenhoven, Rev. Mod. 
Phys. 75, 1 (2003). 

8 H. W. Liu, T. Fujisawa, T. Hayashi, and Y. Hirayama, 
Phys. Rev. B 72, 161305(R) (2005). 

9 V. N. Golovach and D. Loss, Phys. Rev. B 69, 245327 
(2004). 

10 D. T. McClure, L. DiCarlo, Y. Zhang, H.-A. Engel, C. M. 
Marcus, M. P. Hanson, and A. C. Gossard, Phys. Rev. 
Lett. 98, 056801 (2007). 

11 B. Wunsch, M. Braun, J. Konig, and D. Pfannkuche, Phys. 
Rev. B 72, 205319 (2005). 

12 M. R. Graber, W. A. Coish, C. Hoffmann, M. Weiss, J. 
Furer, S. Oberholzer, D. Loss, and C. Schonenberger, Phys. 
Rev. B 74, 075427 (2006). 

13 E. Cota, R. Aguado, and G. Platero, Phys. Rev. Lett. 94, 
107202 (2005). 

14 B. L. Hazelzet, M. R. Wegewijs, T. H. Stoof, and Yu. V. 
Nazarov, Phys. Rev. B 63, 165313 (2001). 

15 J. Fransson, M. Rasander, Phys. Rev. B 73, 205333 (2006). 

16 J. Fransson, Nanotechnology 17, 5344 (2006); New J. 
Phys. 8, 114 (2006). 

17 J. Inarrea, G. Platero and A. H. MacDonald, Phys. Rev. 
B 76, 085329 (2007). 

18 B. Muralidharan and S. Datta, Phys. Rev. B 76, 035432 
(2007). 

19 P. Barthold, F. Hohls, N. Maire, K. Pierz, and R. J. Haug, 
Phys. Rev. Lett. 96, 246804 (2006). 

20 G. Kiefilich, E. Scholl, T. Brandes, F. Hohls, and R. J. 
Haug, Phys. Rev. Lett. 99, 206602 (2006). 

21 J. Aghassi, A. Thielmann, M. H. Hettler, and G. Schon, 
Phys. Rev. B 73, 195323 (2006). 

22 P. Trocha, J. Barnas, Phys. Rev. B 76, 165432 (2007). 

23 I. Weymann, Phys. Rev. B 75, 195339 (2007). 

24 J. N. Pedersen, B. Lassen, A. Wacker, and M. H. Hettler, 
Phys. Rev. B 75, 235314 (2007). 

25 S. A. Wolf, D. D. Awschalom, R. A. Buhrman, J. M. 
Daughton, S. von Molnar, M. L. Roukes, A. Y. Chtchelka, 
and D. M. Treger, Science 294, 1488 (2001). 

26 Semiconductor Spintronics and Quantum Computation, 
ed. by D.D. Awschalom, D. Loss, and N. Samarth 
(Springer, Berlin 2002). 

27 S. Maekawa and T. Shinjo, Spin Dependent Transport in 
Magnetic Nano structures (Taylor & Francis, New York 
2002). 

28 I. Zutic, J. Fabian, S. Das Sarma, Rev. Mod. Phys. 76, 
323 (2004). 



29 S. Maekawa, Concepts in Spin Electronics (Oxford Univer- 
sity Press, 2006). 

30 R. Aguado and D. C. Langreth, Phys. Rev. Lett. 85, 1946 
(2000). 

31 R. Lopez, R. Aguado, and G. Platero, Phys. Rev. Lett. 89, 
136802 (2002). 

32 J. C. Chen, A. M. Chang, and M. R. Melloch, Phys. Rev. 
Lett. 92, 176801 (2004). 

33 D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 
(1998). 

34 D. Loss and E. V. Sukhorukov, Phys. Rev. Lett. 84, 1035 

(2000) . 

35 X. Hu and S. Das Sarma, Phys. Rev. A 61, 062301 (2000). 

36 R. Hanson and G. Burkard, Phys. Rev. Lett. 98, 050502 

(2007) . 

37 M. Julliere, Phys. Lett. A 54, 225 (1975). 

38 J. Barnas and A. Fert, Phys. Rev. Lett. 80, 1058 (1998); 
S. Takahashi and S. Maekawa, Phys. Rev. Lett. 80, 1758 
(1998). 

39 B. R. Bulka, Phys. Rev. B 62, 1186 (2000). 

40 W. Rudzihski and J. Barnas, Phys. Rev. B 64, 085318 

(2001) . 

41 J. Konig and J. Martinek, Phys. Rev. Lett. 90, 166602 

(2003) . 

42 M. Braun, J. Konig, J. Martinek, Phys. Rev. B 70, 195345 

(2004) ; Phys. Rev. B 74, 075328 (2006). 

43 I. Weymann, J. Konig, J. Martinek, J. Barnas, and G. 
Schon, Phys. Rev. B 72, 115334 (2005). 

44 I. Weymann, J. Barnas, J. Konig, J. Martinek, and G. 
Schon, Phys. Rev. B 72, 113301 (2005). 

45 S. Braig and P. W. Brouwer, Phys. Rev. B 71, 195324 

(2005) . 

46 I. Weymann, J. Barnas, Phys. Rev. B 73, 205309 (2006); 
Phys. Rev. B 75, 155308 (2007); Eur. Phys. J. B 46, 289 
(2005); I. Weymann, Europhys. Lett. 76, 1200 (2006). 

47 F. M. Souza, J. C. Egues, and A. P. Jauho, Phys. Rev. B 
75, 165303 (2007). 

48 J. Fransson, Europhys. Lett. 70, 796 (2005). 

49 A. Cottet, W. Belzig, and C. Bruder, Phys. Rev. B 70, 
115315 (2004); Phys. Rev. Lett. 92, 206801 (2004). 

50 A. Cottet and M.-S. Choi, Phys. Rev. B 74, 235316 (2006). 
J. Martinek, Y. Utsumi, H. Imamura, J. Barnas, S. 
Maekawa, J. Konig, and G. Schon, Phys. Rev. Lett. 91, 
127203 (2003); M.-S. Choi, D. Sanchez, and R. Lopez, 
Phys. Rev. Lett. 92, 056601 (2004). 

52 Y. Utsumi, J. Martinek, G. Schon, H. Imamura, and S. 
Maekawa, Phys. Rev. B 71, 245116 (2005). 

53 M. Sindel, L. Borda, J. Martinek, R. Bulla, J. Konig, G. 
Schon, S. Maekawa, and J. von Delft, Phys. Rev. B 76, 
045321 (2007). 

54 I. Weymann and J. Barnas, Phys. Rev. B 77, 075305 

(2008) . 

55 R. Hornberger, S. Koller, G. Begemann, A. Donarini, and 
M. Grifoni, cond-mat/0712.0757 (unpublished). 

56 J. R. Petta, A. C. Johnson, C. M. Marcus, M. P. Hanson, 
and A. C. Gossard, Phys. Rev. Lett. 93, 186802 (2004). 

57 A. C. Johnson, J. R. Petta, and C. M. Marcus, M. P. Han- 
son and A. C. Gossard, Phys. Rev. B 72, 165308 (2005). 

58 H. I. Jorgensen, K. Grove-Rasmussen, J. R. Hauptmann, 
P. E. Lindelof, Appl. Phys. Lett. 89 232113 (2006). 

59 S. Sapmaz, C. Meyer, P. Beliczynski, P. Jarillo-Herrero, 



15 



Leo P. Kouwenhoven, Nano Lett. 6 (7), 1350 (2006). 

60 M. R. Graber, M. Weiss, C. Schonenberger, Semicond. Sci. 
Technol. 21, 64 (2006). 

61 Y. Chye, M. E. White, E. Johnston-Halperin, B. D. Ger- 
ardot, D. D. Awschalom, and P. M. Petroff, Phys. Rev. B 
66, 201301(R) (2002). 

62 H. B. Heersche, Z. de Groot, J. A. Folk, L. P. Kouwen- 
hoven, and H. S. J. van der Zant, A. A. Houck, J. 
Labaziewicz, and I. L. Chuang, Phys. Rev. Lett. 96, 
017205 (2006). 

63 L. Y. Zhang, C. Y. Wang, Y. G. Wei, X. Y. Liu, and D. 
Davidovic, Phys. Rev. B 72, 155445 (2005). 

64 K. Tsukagoshi, B. W. Alphenaar, and H. Ago, Nature 401, 
572 (1999). 

65 B. Zhao, I. M6nch, H. Vinzelberg, T. Miihl, and C. 
M. Schneider, Appl. Phys. Lett. 80, 3144 (2002). 

66 A. Jensen, J. R. Hauptmann, J. Nygard, and P. E. Linde- 
lof, Phys. Rev. B 72, 035419 (2005). 

67 S. Sahoo, T. Kontos, J. Furer, C. Hoffmann, M. Graber, 
A. Cottet, and C. Schonenberger, Nature Physics 1, 102 
(2005). 

68 A. N. Pasupathy, R. C. Bialczak, J. Martinek, J. E. Grose, 
L. A. K. Donev, P. L. McEuen, and D. C. Ralph, Science 
306, 86 (2004). 

69 A. Bernand-Mantel, P. Seneor, N. Lidgi, M. Munoz, V. 
Cros, S. Fusil, K. Bouzehouane, C. Deranlot, A. Vaures, F. 
Petroff, and A. Fert, Appl. Phys. Lett. 89, 062502 (2006). 

70 K. Hamaya, S. Masubuchi, M. Kawamura, T. Machida, M. 
Jung, K. Shibata, K. Hirakawa, T. Taniyama, S. Ishida, 
and Y. Arakawa, Appl. Phys. Lett. 90, 053108 (2007). 

71 K. Hamaya, M. Kitabatake, K. Shibata, M. Jung, M. 
Kawamura, K. Hirakawa, T. Machida, T. Taniyama, S. 
Ishida, and Y. Arakawa, Appl. Phys. Lett. 91, 22107 
(2007); Appl. Phys. Lett. 91, 232105 (2007). 

72 K. Hamaya, M. Kitabatake, K. Shibata, M. Jung, M. 
Kawamura, S. Ishida, T. Taniyama, K. Hirakawa, Y. 



Arakawa, and T. Machida, Phys. Rev. B 77, 081302(R) 
(2008). 

73 H. Yang, S.-H. Yang, S. S. P. Parkin, Nano Lett. 8, 340 
(2008). 

74 B. R. Bulka and T. Kostyrko, Phys. Rev. B 70, 205333 
(2004). 

75 H. Schoeller and G. Schon, Phys. Rev. B 50, 18436 (1994); 
J. Konig, J. Schmid, H. Schoeller, and G. Schon, Phys. 
Rev. B 54, 16820 (1996). 

76 A. Thielmann, M. H. Hettler, J. Konig, and G. Schon, 
Phys. Rev. Lett. 95, 146806 (2005); Phys. Rev. B 68, 
115105 (2003). 

77 D. V. Averin and A. A. Odintsov, Phys. Lett. A 140, 251 
(1989); D. V. Averin and Yu. V. Nazarov, Phys. Rev. Lett. 
65, 2446 (1990). 

78 K. Kang and B. I. Min, Phys. Rev. B 55, 15412 (1997). 

79 D. V. Averin and Yu. V. Nazarov, in Single Charge Tunnel- 
ing, edited by H. Grabert and M. Devoret (Plenum, New 
York 1992). 

80 Ya. M. Blanter and M. Biittiker, Phys. Rep. 336, 1 (2000). 

81 E. V. Sukhorukov, G. Burkard, and D. Loss, Phys. Rev. B 
63, 125315 (2001). 

82 R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, 
L. M. K. Vandersypen, Rev. Mod. Phys. 79, 1217 (2007). 

83 K. Ono and S. Tarucha, Phys. Rev. Lett. 92, 256803 

(2004) . 

84 A. Kogan, S. Amasha, D. Goldhaber-Gordon, G. Granger, 
M. A. Kastner, and H. Shtrikman, Phys. Rev. Lett. 93, 
166602 (2004). 

85 A. W. Holleitner, A. Chudnovskiy, D. Pfannkuche, K. 
Eberl, and R. H. Blick, Phys. Rev. B 70, 075204 (2004). 

86 I. Weymann and J. Barnas, J. Phys.: Condens. Matter 19, 
096208 (2007). 

87 I. Djuric, B. Dong, H. L. Cui, Appl. Phys. Lett. 87, 032105 

(2005) . 



Effects of different geometries on the conductance, shot noise and tunnel 
magnetoresistance of double quantum dots 



OO 

o 
o 

C 

(N 



Ireneusz Weymann 1, 2, Q 

department of Physics, Adam Mickiewicz University, 61-614 Poznan, Poland 

^Department of Theoretical Physics, Institute of Physics, 
Budapest University of Technology and Economics, H-1521 Budapest, Hungary 

(Dated: June 27, 2008) 

The spin-polarized transport through a coherent strongly coupled double quantum dot (DQD) 
system is analyzed theoretically in the sequential and cotunneling regimes. Using the real-time 
diagrammatic technique, we analyze the current, differential conductance, shot noise and tunnel 
magnetoresistance (TMR) as a function of both the bias and gate voltages for double quantum 
dots coupled in series, in parallel as well as for T-shaped systems. For DQDs coupled in series, we 
find a strong dependence of the TMR on the number of electrons occupying the double dot, and 
super-Poissonian shot noise in the Coulomb blockade regime. In addition, for asymmetric DQDs, we 
analyze transport in the Pauli spin blockade regime and explain the existence of the leakage current 
in terms of cotunneling and spin-flip cotunneling-assisted sequential tunneling. For DQDs coupled 
in parallel, we show that the transport characteristics in the weak coupling regime are qualitatively 
similar to those of DQDs coupled in series. On the other hand, in the case of T-shaped quantum dots 
and predict a large super-Poissonian shot noise and TMR enhanced above the Julliere value due to 
increased occupation of the decoupled quantum dot. We also discuss the possibility of determining 
the geometry of the double dot from transport characteristics. Furthermore, where possible, we 
compare our results with existing experimental data on nonmagnetic systems and find qualitative 
agreement. 
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I. INTRODUCTION 

Transport properties of double quantum 
dots (DQDs) have recently attracted much 
interest . 1 - 2 - 3 - 4 - 5 - 6 - 7 - 8 - 9 - 10 - 11 ' 12 - 13 - 14 - 15 ' 16 - 17 - 18 - 19 - 20 - 21 - 22 - 23 - 24 

Since the behavior of double quantum dots resembles 
the behavior of molecules, DQDs are frequently called 
artificial molecules, and are thus considered as ideal 
systems to study the fundamental many-body inter- 
actions between single electrons and spins . 25 i 26 i 27 i 28 i 29 
Double quantum dots exhibit a variety of interesting 
effects, such as for example current rectification due to 
the Pauli spin blockade j^£ i 15 ' 16 i 17 ' 18 negative differential 
conductance' 



pumping fHiii Kondo effect ^-^L^l 



formation of molecular statesj^i 3 .^ S p m 
t) 30 1 3i 1 32 etc Furthermore, 
double quantum dots are also interesting for future ap- 
plications in quantum computingi 33 i 34 i 35 i 36 In addition, 
when the leads are ferromagnetic, transport properties 
of the system strongly depend on the relative orientation 
of the magnetizations of electrodes, leading to the tunnel 
magnetoresistance (TMR) effect, spin accumulation, 
exchange field, eic ^ML^^L^. 

The problem of spin-polarized transport properties has 
been so far mainly addressed in the case of single quan- 
tum dot; . 39 - 40 - 41 - 42 - 43 - 44 - 45 - 46,47 - 48 - 49 - 50,51 - 52 - 53 - 54 In partic- 
ular, in the strong coupling regime it was shown that the 
Kondo peak becomes split when the magnetic configu- 
ration changes from the antiparallel to parallel one, and 
that this splitting can be compensated upon applying ex- 
ternal magnetic fiekh 51 ' 52 i 53 In the weak coupling regime, 
on the other hand, the parity effect of the linear response 



TMR was predicted^ 3 - and the zero-bias anomaly was 
found in the differential conductance when magnetic mo- 
ments of the leads form an antiparallel configuration. 44 
In this paper we extend the existing theoretical studies 
and consider transport through coherent double quantum 
dots weakly coupled to external ferromagnetic leads. We 
note that spin-dependent transport properties of DQDs 
have been analyzed very recently in a few theoretical 
papers i 16 ! 23 ' 55 The considerations were however limited 
to either first-order tunnelin g or the Coulomb blockade 
regime. In particular, in Ref. I16| the formation of a pure 
triplet state was predicted, Ref. [23| deals with transport 
in the deep Coulomb blockade regime, while Ref. [HH 
presents a detailed analysis of sequential transport in the 
case of non-collinearly polarized leads. 

The goal of this paper is thus to analyze the transport 
properties of double quantum dots coupled to collinear 
ferromagnetic leads in the full weak coupling regime, 
i.e. including sequential, cotunneling and cotunneling- 
assisted sequential processes. Furthermore, we also ana- 
lyze the effect of different geometries of the double dots 
on transport, in particular, the cases of DQDs coupled in 
series, in parallel, as well as T-shaped double quantum 
dots are considered. The comparison of numerical results 
obtained for different geometries may in principle help in 
determining the system's geometry, which may be of im- 
portance in discussing and understanding experimental 
results, especially on self-assembled quantum dots. The 
present analysis is based on the real-time diagrammatic 
technique which consists in a perturbation expansion of 
the density matrix of the system under consideration, 



and the relevant operators, with respect to the coupling 
to external leads. The advantage of using the real-time 
diagrammatic technique is that it takes into account the 
effects of the exchange field and the renormalization of 
the dot levels in a fully systematic way. In particular, 
in the case of asymmetric DQDs coupled in series, we 
consider transport in the Pauli spin blockade regime and 
show that the leakage current in the blockade region re- 
sults from cotunneling and cotunneling-assisted sequen- 
tial tunneling processes. As far as the shot noise is con- 
cerned, we show that the noise is super-Poissonian in the 
Coulomb blockade regime and drops to sub-Poissonian 
value in the sequential tunneling regime, where the Fano 
factor approaches one half. This behavior is observed 
for parallel and serial DQD's geometries. In turn, for 
T-shaped double quantum dots, we find a large super- 
Poissonian shot noise due to increased occupation of the 
decoupled dot. We also analyze the TMR effect and find 
its strong dependence on the transport regime and num- 
ber of electrons in the DQD in the ground state. For 
DQDs coupled in series and in parallel the TMR takes 
the values ranging from around a half of the Julliere TMR 
to its full valuer 7 , whereas for T-shaped DQDs the TMR 
may be enhanced above the Julliere value. 

The systems considered in this paper may be realized 
experimentally for example in lateral and vertical semi- 
conductor quantum dot s 6 ' 8 ' 10 ! 56 ' 57 or single wall metal- 
lic carbon nanotubes with top gate electrodes i 12 i 58 ' 59 i 60 
The latter systems are of particular interest because, by 
tuning the gates, it is possible to change the charge 
on each dot separately. Unfortunately, the aforemen- 
tioned experiments concerned only DQDs coupled to 
nonmagnetic leads. There are several experimental real- 
izations of single quantum dots attached to ferromagnetic 
le^^LS^S^M^SMmiJMS. while experimental 
data on spin-polarized transport through double quan- 
tum dots are lacking. We thus believe that the results 
presented in this paper will be of assistance in discussing 
future experiments. 

The paper is organized as follows: The model of the 
considered system is presented in Sec. II, whereas the 
method employed in calculations is described in Sec. III. 
In Section IV we present and discuss the numerical results 
for double quantum dots coupled in series. In this section 
we also focus on the role of second-order processes in 
transport. In addition, we also analyze transport in the 
Pauli spin blockade regime. In Sec. V we deal with 
DQDs coupled in parallel, whereas in Sec. VI we consider 
transport through T-shaped quantum dots. Finally, the 
conclusions are given in Sec. VII. 



II. MODEL 

The schematic of a double quantum dot coupled to fer- 
romagnetic leads is shown in Fig. [TJ It is assumed that 
the magnetizations of the leads are oriented collinearly, 
so that the system can be either in the parallel or an- 
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FIG. 1: (Color online) Schematic of a double quantum dot 
weakly coupled to external ferromagnetic electrodes. The 
magnetizations of the leads can form either parallel or an- 
tiparallel magnetic configuration. The first (j — 1) and sec- 
ond (j = 2) dots are coupled to each other via the hopping t, 
and to the left (r — L) and right (r = R) leads with the cou- 
pling strength T r j . The double quantum dot is assumed to be 
symmetrically biased. By adjusting the couplings, the system 
can smoothly cross over from the serial to parallel geometry. 

tiparallel magnetic configuration. The Hamiltonian of 
the system is given by 

H = Hl + + #DQD + Ht (1) 

The first two terms describe noninteracting itinerant elec- 
trons in the leads, H r = J2ka e rkacl ka c r ] icr for the left 
(r = L) and right (r = R) lead, where e r ko- is the energy 
of an electron with the wave vector k and spin a in the 
lead r, and c^ k(T (c r kcr) denotes the respective creation 
(annihilation) operator. The double dot is described by 
the Hamiltonian 

j=l,2 a j=l,2 

+ U'J2 ni„n 2a > + t ^2(4 a d 2a + 4 a d la ) ,(2) 

with rij a — d^dju, where (dj a ) is the creation (anni- 
hilation) operator of an electron with spin a in the first 
(j = 1) or second (j = 2) quantum dot, and e 3 - is the 
corresponding single-particle energy. The Coulomb in- 
teraction on the first and second dot is assumed to be 
equal and is described by U, while U' corresponds to the 
inter-dot Coulomb correlation. The last term of Hdqb 
describes the hopping between the two dots with t being 
the hopping parameter. We assume that the hopping pa- 
rameter is large, so that there is a considerable overlap 
of the wave functions of the two dots, leading to the for- 
mation of molecular many-body states, through which 
transport takes place. In addition, we also note that 
an exchange interaction between spins in the two dots 
may lead to the formation of singlet and triplet states*^ 
However, experimentally, this exchange interaction was 
found to be rather small as compared to the other en- 
ergy scales^ therefore, in the following considerations, 
we will neglect it. 

The tunneling processes between the DQD and elec- 
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trodes are described by the Hamiltonian, 

r=L,Rj=l,2 k<T 

where t r j denotes the tunnel matrix elements between 
the rth lead and the jth dot. The coupling of the jth 
dot to the rth lead can be written as T°j = 2n\t r j\ 2 p^, 
where p% is the spin-dependent density of states of lead 
r. By introducing the definition of the spin polarization 
of lead r, p r = (p+ — p~)/(p+ + Pr)i one can write, 

r 5 _) = r rj (i±p r ), with r rj = (r+.+r-.)/2, where r+ 

and r~- describe the coupling of the jth dot to the spin- 
majority and spin-minority electron bands of the lead r, 
respectively. 

Because the double dot Hamiltonian, Eq. is not 
diagonal in the local basis, we perform a unitary trans- 
formation, W HdqdII = Hbqd, to a new basis in which 

#dqd is diagonal, #dqd|x) = £ x \x)- Tne eigenvectors 
\x) are the many-body states of the double quantum 
dot, while the eigenvalues e x denote the corresponding 



energies 
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III. METHOD 

In order to calculate the spin-dependent transport 
properties of a double quantum dot in the sequential 
and cotunneling regimes, we employ the real-time dia- 
grammatic techniquei 43 ' 75 i 76 It consists in a systematic 
perturbation expansion of the reduced density matrix of 
the considered system and the current operator with re- 
spect to the dot-lead coupling strength T. The time evo- 
lution of the reduced density matrix can be visualized 
as a sequence of irreducible self-energy blocks, W xx > , on 
the Keldysh contour, where W xx ' describe transitions be- 
tween the many-body DQD states \\) and \x')- The ele- 
ments W xx > set up a self-energy matrix W. Within the 
matrix notation introduced by Thielmann et al.j^ the 
stationary elements of the reduced density matrix can be 
found from 



(Wp 



si \ 

h 



r<5. 



(4) 



where p st is the vector containing probabilities and the 
matrix W is the modified matrix W so as to include 
the normalization of probabilities. The current flowing 
through the system can be then found from 



2fi 



Tr{WV } 



(5) 



The matrix W 1 is the self-energy matrix with one inter- 
nal vertex resulting from the expansion of the tunneling 
Hamiltonian replaced by the current operator. 

To calculate the transport properties order by order in 
tunneling processes, we expand the self-energy matrices, 
W = W« + W (2 ) + . . . , W 1 = W 1 ^) + W 1 ' 2 ) + . . . , and 



the dot occupations, p st = p st (°) + p st (!) + . . . , respec- 
tively. Then, the first-order (sequential) and the second- 
order (cotunneling) currents are given by 



j(D 



2H 

e 



Tr { W I(l) p st(0) } 



jP) = ^{w^p-M+W^yW}. (6) 

On the other hand, the zeroth and first-order occupation 
probabilities can be found from the following equations 



r<s v 



j x — _L UxXo , 

w<V t(1) +w(V t (°> - 0, (7) 

where W« [w( 2 )j is given by W« [W^ 2 )] with one 

arbitrary row \o replaced by (r, . . . , T) [(0, . . . , 0)] due 

to the normalization of probabilities, Yl x Px = ^™' ' 

As one can see from the above formulas, to deter- 
mine the transport properties it is necessary to calcu- 
late the elements of the corresponding self-energy ma- 
trices. This can be done using the respective diagram- 
matic rules i 43 i 75 ' 76 Although the first-order calculation 
is rather simple, the general analytical formulas for the 
self-energies in the second order are rather complicated 
due to many virtual states through which the cotunneling 
processes can take place. In the Appendix A, as an exam- 

(2) 

pie, we present the contribution coming from I^vu WiV) 
where N is the charge state of the double dot. 

In addition, in the present paper we also analyze the 
zero-frequency current noise of the double quantum dot 
system, S = 2j^ dfc((/(t)/(0) + J(0)/(t)) - 2(/) 2 ). For 
low bias voltages the current noise is dominated by ther- 
mal noise, while for \eV\ > k-&T the noise associated with 
the discrete nature of charge (shot noise) dominates. 80 
Within the real-time diagrammatic technique, a general 
formula for the shot noise has been derived in Ref. [76| 
taking into account the non-Markovian effects. The shot 
noise in the first-order can be found from the correspond- 
ing expression 



- ^Tr 

h 



(\y II(1) + W I(1) P (-1) W I(1) ) p st(0) 



while the cotunneling current noise is given by 



(8) 



5(2) = £_Tr(|w II(2) + W I(2) P ( ~ 1} W I(1) 

+ w I(1) P (_1) W I(2) + W I(1) P (0) W I(1) 
+ W^Q^AWWl p st (°) 

+ [w^ + W'WpHlwiWlp^D^ (9) 

where QW = p st (") ® e T , with e T = (1, . . . , 1). The ob- 
jects P ( ~ 1} and P (0) are given by, W^pt" 1 ) = QW - 1, 
and, WWpM +W( 2 )P(-!) = 1(QW -aW^QW), re- 
spectively, with 1 being the unit vector with row xo set to 
zero. On the other hand, the matrices W 11 ^ 1 ) and W 11 ' 2 ' 



4 




are the first and second-order self-energy matrices with 
two internal vertices replaced by the current operator, 
while <9W« and dW 1 ^ are partial derivatives of W' ' 
and W ! W with respect to the convergence factor of the 
Laplace transform^ 

By taking into account all the first-order and second- 
order contributions to the self-energies, we are able to 
resolve the transport properties in the full range of the 
bias and gate voltages. The first order of expansion with 
respect to the coupling corresponds to sequential tunnel- 
ing, while the second order is associated with cotunncl- 
ing. Sequential tunneling is allowed if the applied bias is 
larger than the threshold voltage, i.e. when the energy 
provided by the transport voltage is comparable with the 
charging energy. Otherwise, the system is in the Coulomb 
blockade regime where sequential tunneling is exponen- 
tially suppressed and the current flows due to cotunncl- 
ing, which involves correlated tunneling through virtual 
states of the system i 77 ' 78 ' 79 

Among different cotunncling processes one can gen- 
erally distinguish two types of processes: non-spin-flip 
ones, which do not affect the DQD state, and the spin- 
flip ones. Furthermore, one can also have double-barrier 
cotunncling, which contributes directly to the current, 
and single-barrier cotunneling, which affects the double 
dot occupations and, thus, indirectly, the current. 



IV. DOUBLE QUANTUM DOTS COUPLED IN 
SERIES 

In this section we present and discuss numerical re- 
sults on double quantum dots coupled in scries. In order 
to realize this geometry we set IY2 = Tri = and as- 
sume Tli = Tri2 = r/2, see Fig. [T] Furthermore, in 
the following we will also distinguish between symmetric 
and asymmetric DQDs. In the former case, the level po- 
sition of each dot is the same, E\ = £2, while in the latter 
case the levels are detuned, £\ ^ £2- In addition, in this 
section we will emphasize the role of second-order pro- 
cesses in transport. In order to ascribe observed features 
to respective tunneling processes, we will therefore also 
present results obtained within the sequential tunneling 
approximation, i.e. when considering only the first-order 
tunneling processes. 



A. Symmetric double quantum dots 

The differential conductance as a function of the bias 
voltage and the position of the dot levels for the parallel 
and antiparallel configurations in shown in Fig. Ufa) and 
(b). Because experimentally the position of the dots' lev- 
els can be changed upon applying a gate voltage, Fig. [3] 
effectively shows the bias and gate voltage dependence of 
the conductance. First of all, it can be seen that the dif- 
ferential conductance displays a characteristic Coulomb 
diamond structure. The diamonds in the low bias volt- 
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FIG. 2: (Color online) The differential conductance in the par- 
allel Gp (a) and antiparallel Gap (b) magnetic configurations 
as a function of the bias voltage V and the position of the dots' 
levels e = ei — E2 for double quantum dots coupled in series. 
The parameters are: fcp,T = 0.15 meV, U — 2 meV, U' = 1 
meV, t = 0.25 meV, r L2 = T R i = 0, lYi = r R2 = T/2, with 
r = 0.1 meV, and p = 0.5. Because the double quantum dot 
is symmetric, e\ = £2, there is no asymmetry associated with 
the bias reversal. For comparison in part (c) we also show 
the density plot of the differential conductance in the parallel 
configuration calculated using only the first-order tunneling 
processes, Gp . 



age regime correspond to the Coulomb blockade regime 
where the sequential current is exponentially suppressed 
and the current flows due to cotunneling. With increas- 
ing the bias voltage the excited states start participating 
in transport, which leads to additional lines in the dif- 
ferential conductance. The features described above are 
rather associated with the energy spectrum and charge 
states of the DQD than with the ferromagnetism of the 
leads, therefore they are present in both magnetic con- 
figurations. The spin dependence stemming from ferro- 
magnetic electrodes gives rise to a difference of the con- 
ductance in the parallel and antiparallel configurations, 
see Fig. [^a) and (b). The conductance in the parallel 
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FIG. 3: (Color online) The total (first and second order) TMR 
(a) and the first-order TMR (b) as a function of the bias 
voltage V and the level position e for parameters the same as 
in Fig. [2] The two figures are plotted in the same scale. 



configuration is generally larger than that in the antipar- 
allel configuration. This is due to the spin asymmetry 
in the couplings when the leads' magnetizations are an- 
tiparallel. As in the parallel configuration the majority 
(minority) electrons of, let us say, left lead tunnel to the 
majority (minority) electron band of the right lead, in 
the antiparallel configuration the situation is reversed - 
the majority (minority) electrons of the left lead tun- 
nel to the minority (majority) electron band of the right 
lead. This leads to the difference of the conductance in 
both magnetic configurations, see Fig. [2j and to the cor- 
responding tunnel magnetoresistance effect. In addition, 
in Fig. [^c) we also present the bias and gate voltage de- 
pendence of the differential conductance for the parallel 
configuration calculated within the sequential tunneling 
approximation, Gp\ By comparing Fig. dja) and (c), 
one can see that cotunneling leads to broadening of the 
resonance peaks, lowering the magnitude of the differen- 
tial conductance, and gives rise to finite conductance in 
the Coulomb blockade regimes. 

The bias and gate voltage dependence of the TMR is 
shown in Fig.[3ja), while Fig.[3jb) presents the first-order 
TMR, TMR (1) . Intuitively, the TMR is a measure of the 
system's transport properties change when the magnetic 
configuration switches from the parallel to the antiparal- 
lel one. It is defined as2L2&^2 



TMR = 



/At 



(10) 




where Ip (Zap) is the current flowing through the sys- 



FIG. 4: (Color online) The linear conductance in the parallel 
(solid line) and antiparallel (dashed line) magnetic configura- 
tions (a) and the linear TMR (b) as a function of the level 
position e for parameters the same as in Fig. [2] For com- 
parison, the dotted curves show the linear conductance and 
the TMR calculated taking into account only the first-order 
processes. 



tern in the parallel (antiparallel) configuration. By com- 
parison with Fig. [21 one can easily identify the different 
transport regimes. Furthermore, by comparing the total 
and sequential TMR, one can immediately see that co- 
tunneling modifies the TMR mainly in the low bias volt- 
age regime, i.e. in the Coulomb blockade regime and in 
the regime where sequential tunneling is suppressed due 
to the absence of levels in the energy window provided 
by transport voltage. In order to discuss and see more 
clearly the behavior of the TMR on applied voltages, in 
Fig. 0] we show the linear response conductance and re- 
sulting TMR, while in Fig. [5] we display the transport 
properties in the nonlinear response regime. 

The linear conductance in both magnetic configura- 
tions is shown in Fig. 0Ja). When lowering the position 
of the DQD levels the conductance displays four reso- 
nance peaks associated with subsequent occupation of 
the corresponding charge states. The dotted curves in 
Fig. 2] present results obtained in the first-order approx- 
imation. It can be seen that off resonance the current is 
mainly dominated by cotunneling. On the other hand, 
on resonance the sequential processes dominate current, 
while cotunneling only slightly affects the conductance, 
leading to renormalization of the DQD levels and, thus, 
slightly shifting the position of the conductance peaks, 
see Fig. Ufa). Interestingly, the second-order processes 
have a rather large impact on the linear TMR shown in 
Fig. IUb). First of all, the linear response TMR exhibits 
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a strong dependence on the gate voltage, i.e. on the 
number of electrons in the double quantum dot. When 
the DQD is either empty or fully occupied, the TMR 
is given by the Julliere value, TMR = 2p 2 /(I — p 2 )£L 
This is due to the fact that in those transport regimes 
the current is driven by elastic cotunneling processes 
which do not affect the DQD state in any way. Thus, 
as far as the TMR is concerned, the system behaves as a 
single ferromagnetic tunnel junction^ However, in the 
Coulomb blockade regions, the TMR becomes generally 
suppressed due to the presence of inelastic cotunneling 
which introduces spin-flip processes in the system. The 
dotted line in Fig.fJJb) shows the TMR calculated within 
the sequential tunneling approximation - it is given by 
TMR = p 2 /(l - p 2 )M& At this point, we would like 
to note that the results for the linear conductance cal- 
culated by including only the first-order processes are 
rather reliable in the whole range of e. On the other 
hand, the results for the linear response TMR are com- 
parable to those obtained within the sequential tunneling 
approximation only on resonance, where sequential pro- 
cesses dominate, while off resonance they are completely 
unreliable, see Fig. [¥] Therefore, in order to properly 
analyze the dependence of the TMR in the full range 
of parameters one has to take into account cotunneling 
processes. 

The bias voltage dependence of the current, differen- 
tial conductance, tunnel magnetoresistance and the Fano 
factor is shown in Fig. [5] for e = —2 meV and e = 2 
meV. The first situation, e — — 2 meV, corresponds to 
the case when the ground state of the double quantum 
dot is doubly occupied, while the second one, e — 2 meV, 
corresponds to the case of empty double dot. In all cases, 
due to the spin asymmetry of tunneling processes, the 
current in the parallel configuration is larger than the 
current in the antiparallel configuration. In addition, the 
i-v curves display characteristic Coulomb steps. With 
increasing the bias voltage more and more charge states 
start to participate in transport which gives rise to the 
corresponding steps in the current-voltage characteris- 
tics, and peaks in the differential conductance. Further- 
more, the single-electron charging effects also lead to the 
oscillatory-like behavior of the TMR effect, see Fig. [5jc) 
and (g). In the Coulomb blockade regime, e = — 2 meV, 
the TMR is much suppressed as compared to the case of 
e = 2 meV. This is due to the presence of spin-flip cotun- 
neling processes in the case of the doubly occupied DQD, 
as already discussed in association with Fig. |U However, 
when increasing the bias voltage the TMR increases and, 
at the threshold for sequential tunneling, reaches a local 
maximum, see Fig.[5][c). This effect is due to the nonequi- 
librium spin accumulation which is induced in the DQD 
system with increasing the transport voltage. (Similar 
effect can be observed when the double dot is singly oc- 
cupied, i.e. when e — —0.5 meV.) On the other hand, 
when the DQD is empty, the TMR at low bias voltage 
acquires the Julliere valued and then, with increasing 
V, becomes generally suppressed, see Fig. [SJg). 




FIG. 5: (Color online) The current in units of Io ~ eT/h (a,e), 
differential conductance (b,f) and the Fano factor (d,h) in 
the parallel (solid line) and antiparallel configuration (dashed 
line) as well as the TMR (c,g) as a function of the bias voltage 
for different values of the level position e as indicated in the 
figure. £ = — 2 meV (a-d) corresponds to the case when the 
ground state of the DQD is doubly occupied, while for e = 2 
meV (e-h) the double dot is empty. The parameters are the 
same as in Fig. [2] The dotted curves show the results obtained 
in the sequential tunneling approximation. 



In addition, in Fig. O we also show the Fano factor, 
F = Sy(2e|I|), calculated in both magnetic configura- 
tions and for two values of e, as indicated in the fig- 
ure. First of all, one can see that the Fano factor be- 
comes divergent at low bias voltages. This is due to 
the thermal noise which dominates the current noise as 
V — > 0, while 1^0, leading to F — > cc<££ Further- 
more, it can be seen that the Fano factor in the Coulomb 
blockade regime is slightly larger in the parallel config- 
uration than in the antiparallel one. This is associated 
with an additional positive cross-correlations, which con- 
tribute to the current noise, due to the ferromagnetism 
of the electrodes. In the Coulomb blockade regime, there 
is a larger asymmetry between the cotunneling processes 
through the majority and minority spin channels in the 
parallel configuration than in the antiparallel one. This 
effectively increases the current fluctuations in the par- 
allel configuration and leads to the corresponding dif- 
ference in the Fano factors, see Fig. E^d). On the other 
hand, in the sequential tunneling regime, the information 
about the spin asymmetry of tunneling processes between 
the DQD and the corresponding lead is rather contained 
in the nonequilibrium spin accumulation induced in the 
double dot. This spin accumulation is larger in the an- 
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FIG. 6: (Color online) The total Fano factor in the parallel 
(a) and antiparallel (b) magnetic configurations as a function 
of the bias voltage V and the position of the dots' levels e 
for parameters the same as in Fig. [2] Part (c) presents the 
Fano factor in the parallel configuration calculated in the first- 
order approximation. Because the Fano factor is divergent 
for \eV\ < kuT, this transport regime is denoted with a thick 
black line. 



tiparallel configuration, yielding generally an enhanced 
Fano factor in the antiparallel configuration as compared 
to the parallel one. Moreover, we also note that in the 
Coulomb blockade regime, e = —2 meV, the shot noise 
becomes super-Poissonian (F > 1) and drops to the sub- 
Poissonian value (F < 1) at the threshold for sequential 
tunneling, see Fig.[5jd). This super-Poissonian shot noise 
is associated with bunching of inelastic cotunncling pro- 
cesses through the system.! 54 ' 81 On the other hand, in 
the case of empty DQD, e = 2 mcV, in the cotunnel- 
ing regime the current flows due to elastic second-order 
processes which obey the Poissonian statistics. In this 
case the shot noise is given by S = 2e\I\ and the Fano 
factor is simply equal to unity, see Fig.[5Jh), irrespective 



of magnetic configuration of the system. As concerns 
the sequential tunneling regime, the shot noise is then 
generally sub-Poissonian, irrespective of the position of 
the DQD level e. This indicates the role of correlations 
in electronic transport, in particular, the Coulomb cor- 
relations and charge conservation^ Our results are in 
qualitative agreement with experimental data on tunnel- 
ing through vertically coupled self-assembled InAs quan- 
tum dots where also super-Poissonian shot noise has been 
founds 

Finally, to make the present analysis self-contained, we 
also show the bias and gate voltage dependence of the to- 
tal Fano factor for the parallel and antiparallel magnetic 
configuration, see Fig.^a) and (b), respectively. In addi- 
tion, for comparison, the Fano factor calculated by taking 
into account only first-order tunneling processes is de- 
picted in Fig. [6](c). The black lines around the zero bias, 
\eV\ < ksT, mark the transport regime where the Fano 
factor is divergent due to finite thermal noise. The differ- 
ent behavior of the shot noise is now clearly visible. In the 
Coulomb blockade regime the noise is super-Poissonian, 
in the cotunneling regime when the DQD is either empty 
of fully occupied (current is mediated by elastic cotun- 
neling) the shot noise is Poissonian, and in the sequential 
tunneling regime the noise drops to sub-Poissonian value. 



B. Asymmetric double quantum dots: Pauli spin 
blockade 

By applying a gate voltage to each quantum dot, it 
is possible to tune the dot levels separately. So far, we 
have considered the case of symmetric DQD, i.e. when 
£\ = £2- In this situation the transport characteristics 
were symmetric with respect to the bias reversal. This 
is however not the case for s\ =/= £2, where the current 
becomes asymmetric with respect to the bias reversal, 
leading to the Pauli spin blockade and negative differen- 
tial conductance, as observed experimentally 6,8 In the 
following, assuming realistic parameters of the double 
quantum dot system, we analyze transport properties in 
the regime where the Pauli spin blockade effects are vis- 
ible. 

The mechanism leading to the spin blockade was theo- 
retically discussed by Fransson et al^ and Muralidharan 
et al^ These considerations were however restricted only 
to the first-order tunneling processes, which dominate the 
current out of the blockade regime. As shown experimen- 
tally by Ono et al.^ in the spin blockade regime there is a 
finite leakage current, which cannot be explained within 
the sequential tunneling approximation. To explain the 
existence of this leakage, it was proposed that non-zero 
current in the Pauli blockade is associated with spin-flip 
processes induced by hyperfine interaction^ In the fol- 
lowing, we show that the leakage current results just from 
the interplay between different intrinsic tunneling pro- 
cesses driving the current. The hyperfine interaction, or 
coupling to a phonon bath, may, of course, increase the 



8 




-3-2-10123 
V(mV) 



FIG. 7: (Color online) The current (a) and Fano factor (c) in 
the parallel (solid) and antiparallel (dashed) magnetic config- 
urations as well as the TMR (b) as a function of the bias 
voltage for DQD coupled in series. The parameters are: 
k B T = 0.05 meV, ei = -1 meV, e 2 = -2 meV, A = 3 meV, 
U = 2 meV, U' = 1 meV, t = 0.05 meV, where A is the level 
spacing in each dot. The other parameters are the same as in 
Fig. [2] The inset in part (a) displays the current in the paral- 
lel magnetic configuration in the Pauli spin blockade regime. 
The dotted lines correspond to the first-order calculation. 



leakage, but is not necessary for the observation of a finite 
current in the Pauli spin blockade regime. 

When the DQDs levels are detuned, £2 < £1, and once 
the first and second dot become singly occupied, the cur- 
rent may be suppressed in some range of the bias voltage. 
This is associated with the full occupation of two-electron 
triplet states of the DQD^ The Pauli spin blockade may 
be lifted when the applied bias voltage admits another 
(exited) charge states to participate in transport. On the 
other hand, when the voltage is reversed, it is energeti- 
cally allowed that the electron from the first dot tunnels 
to the left lead and then another electron from the right 
lead enters the DQD, lifting the Pauli spin blockade. 



This can be seen in Fig. [7] where we present the bias 
voltage dependence of the current, Fano factor and the 
TMR effect. In order to make the calculations more re- 
alistic and to allow for excited states of the system, we 
have now taken into account four different orbital levels, 
two in each dot. Such system is described by straight- 
forward extension of the hamiltonian, Eq. where the 
level spacing in the dots is described by the parameter A. 
First of all, one can see that the transport characteristics 
are asymmetric with respect to the bias reversal. This is 
associated with the fact that by detuning the DQDs lev- 
els the symmetry of transition rates between the left and 
right leads has been broken. Furthermore, for positive 
bias voltage the current is suppressed in a broad range of 
the bias voltage (0 < eV < U) due to the full occupation 
of the triplet states, while for negative bias the Pauli spin 
blockade is lifted. In addition, the blockade can also be 
lifted when the positive bias voltage is increased further, 
eV > U, so that tunneling through the second level of 
the second dot is allowed, see Fig. [TJa). 

In the inset of Fig. EJa) we show the current in the 
parallel configuration just in the spin blockade regime. 
For comparison, we also display the current calculated 
within the sequential tunneling approximation. One can 
see that the sequential current is suppressed as compared 
to the total current (calculated taking into account co- 
tunneling processes). The suppression of the sequential 
current is generally governed by the ratio 2t/|e 2 — ei\M> 
When 2t/\s2— £i| <C 1, the occupation of the two-electron 
triplet states approaches unity and the first-order current 
becomes fully blocked. However, the second-order pro- 
cesses are still allowed, leading to a finite leakage current 
in the spin blockade regime. One can distinguish different 
contributions coming from cotunneling: 

(i) The double-barrier elastic second-order processes 
which contribute directly to the current. 

(ii) The double-barrier spin-flip cotunneling which con- 
tributes to the current and reduces the occupation of the 
triplet state, this way opening the DQD for the sequential 
tunneling. 

(iii) The single-barrier spin-flip cotunneling which does 
not contribute directly to the current but flips the spin 
in the DQD and allows the sequential processes to occur. 

As a consequence, in the Pauli spin blockade regime 
the current flows due to cotunneling and spin-flip 
cotunncling-assisted sequential tunneling. The key role 
is played by the spin-flip processes which decrease the 
occupation of the two-electron triplet states and lead to 
a finite occupation of singlet states. The importance of 
the spin-flip processes in explanation of the finite leakage 
current has been invoked in Ref. [171 ] . The leakage was 
there associated with spin-flip processes induced by hy- 
pcrfine interaction through the Overhauser effect. Here, 
we show that the leakage results just from the interplay 
between different tunneling processes, i.e. just from the 
pure nature of tunneling processes. We also note that the 
effects of cotunneling are rather independent of the mate- 
rial from which the dots are built, unlike the Overhauser 
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FIG. 8: The total (sequential plus cotunneling) current (solid 
line) and the sequential current (dotted line) in the paral- 
lel configuration as a function of the hopping between the 
two dots calculated in the middle of the Pauli spin blockade 
regime, V — 1 mV. The parameters are the same as in Fig. [7] 
The dependence of the current on 2t/\e2 — ei| in the antipar- 
allel configuration is qualitatively similar. 



fields Thus, for example, in GaAs DQDs the leakage 
current in the Pauli spin blockade would be associated 
with both cotunneling and hyperfine contributions, while 
in carbon nanotube DQDs it would be mainly due to co- 
tunneling. Furthermore, as shown in experiments by Ono 
and Tarucha, 83 who studied the dependence of the leak- 
age current on the applied magnetic field, the hyperfine 
interaction starts to play a role at certain finite in-plane 
magnetic fields, at which the nuclei become polarized, 
leading to a sudden jump of the leakage current as a 
function of magnetic field. This suggests that the main 
contribution to the current in the spin blockade regime 
in the absence of magnetic field may come from cotun- 
neling. In the weak coupling regime, typical values of the 
dot-lead coupling strength T are of the order of /ieV, 84 
which gives the leakage current of the order of 10 _3 Jo « 1 
pA (see Fig. [8]), in agreement with experimental results^ 3 - 
Finally, we notice that to explain the jump in the leakage 
current when sweeping magnetic field, besides cotunnel- 
ing, one also has to include the hyperfine interaction. 
Nevertheless, this goes beyond the scope of the present 
paper. 

In Fig. [8] we present the logarithmic dependence of the 
sequential and total currents in the parallel configura- 
tion on the hopping parameter t. It is clearly visible that 
for 2i/|e 2 — £i | <C 1, the sequential current is smaller 
by two orders of magnitude than the total current. On 
the other hand, when 2t/\e2 — £i\ ~ 1, the sequential 
current becomes of the same order as the total current, 
but still, noticeably, the total current is larger than that 
calculated using only the first-order tunneling processes. 
The saturation of the two currents for 2i/|e 2 — £i| -C 1, 
see Fig. HI is associated with a finite temperature. How- 
ever, because sequential tunneling depends exponentially 
on temperature while cotunneling only algebraically, at 
lower T the difference between the two currents would be 



even more pronounced. The dependence of the current 
in the antiparallel configuration is qualitatively similar to 
the one shown in Fig. [8l We also note that generally the 
above mentioned mechanism responsible for the leakage 
current does not depend on whether the leads are made 
of ferromagnetic or nonmagnetic material. 

The TMR as a function of the bias voltage is displayed 
in Fig. [Tfb). One can see that the TMR calculated us- 
ing the first-order processes and that calculated taking 
into account cotunneling are generally similar, except for 
the Pauli spin blockade regime. In this transport regime 
the sequential TMR is much overestimated. To under- 
stand this behavior we note again that the blockade re- 
sults from the full occupation of DQD's triplet states. 
The probability is equally distributed between the three 
components of the triplet in the parallel configuration. 
However, in the antiparallel configuration, it turns out 
that the charge is mainly accumulated in the S z = 1 
component of the triplet. As a consequence, in the par- 
allel configuration all the three components of the triplet 
participate in transport, while in the antiparallel config- 
uration only one. Because the sequential current in the 
blockade regime is mainly associated with thermal fluc- 
tuations, the difference in the number of states relevant 
for transport leads to large sequential TMR in blockade 
regime. This difference also gives rise to the correspond- 
ing difference between Fano factors calculated in the se- 
quential tunneling approximation, see the dotted curves 
in Fig.[7Jc). The shot noise is then super-Poissonian and 
larger in the parallel configuration. However, as pointed 
above, the leakage current in the Pauli spin blockade 
regime is due to cotunneling and cotunneling-assisted 
sequential tunneling processes. The total current flows 
not only due to thermal fluctuations but due to corre- 
lated tunneling through virtual states of the DQD. This 
fact generally decreases the current fluctuations and the 
difference between the two magnetic configurations. As 
a result, the total TMR and Fano factor become sup- 
pressed as compared to the sequential tunneling results, 
see Fig. [7Jb) and (c), respectively. In addition, we also 
notice that, irrespective of the magnetic configuration of 
the system, the shot noise is super-Poissonian in the spin 
blockade regime and drops to sub-Poissonian value out of 
the blockade regime, see Fig.JT^c). The super-Poissonian 
shot noise is generally due to bunching of spin-flip co- 
tunneling and cotunneling-assisted sequential tunneling 
processes. 



V. DOUBLE QUANTUM DOTS COUPLED IN 
PARALLEL 

In this section we analyze the spin-polarized transport 
through double quantum dots coupled in parallel. This 
geometry can be realized by setting T T j = T/2, for r — 
L, R and j = 1, 2, see Fig. [T] We note that in the case of 
negligible hopping between the two quantum dots, t — > 
0, the behavior of the DQD system resembles that of a 
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single multi-level quantum dot. The problem of spin- 
dependent transport through multi-level quantum dots 
has been addressed very recently in Ref . [5J] and will not 
be considered here. 

In Fig. [9] we present the density plots of the differential 
conductance and Fano factor in the parallel configuration 
as well as of the TMR effect as a function of the bias 
and gate voltages. The behavior of the conductance and 
the shot noise in the antiparallel configuration is qualita- 
tively similar as in the parallel configuration, therefore, it 
is not shown here. The information about the difference 
in transport in the two magnetic alignments is contained 
in the TMR, whose magnitude reflects the asymmetry 
in tunneling processes when the leads are parallel or an- 
tiparallel to each other. 

The differential conductance, shown in Fig. [Ha), dis- 
plays characteristic Coulomb diamonds. Because the in- 
ternal energy structure of the DQD is generally the same 
as in the case of DQD coupled in series, the differential 
conductances are qualitatively similar, see Figs. [2a) and 
UK a). The main difference is that in the case of parallel 
DQDs the level of each dot is coupled both to the left 
and right leads, which generally results in an enhanced 
conductance for parallel DQDs as compared to DQDs 
coupled in series. Similarly the TMR shown in Fig.[9jb), 
its bias and gate voltage dependence is only slightly mod- 
ified as compared to Fig. [3ja). Again, the linear TMR 
exhibits a strong dependence on the occupation number 
of the DQD, with TMR given by the Julliere formula 
for empty and fully occupied double quantum dot, and 
is much suppressed in other transport regimes due to 
spin-flip cotunneling. In the Coulomb blockade regime, 
the transport properties are mainly conditioned by the 
spin-flip cotunneling, which gives a dominant contribu- 
tion to the current. Furthermore, the bunching of spin- 
flip cotunneling in the Coulomb blockade regime leads to 
super-Poissonian shot noise, see Fig. [5Jc). The super- 
Poissonian shot noise is however more pronounced in 
Coulomb blockade regimes with an odd number of elec- 
trons in the double dot. On the other hand, for voltages 
above the threshold for sequential tunneling, the shot 
noise becomes sub-Poissonian, with the Fano factor ap- 
proaching 1/2. 

Summing up, we notice that generally in the weak 
coupling regime the effect of different geometries of the 
double dot system leads only to qualitative difference in 
transport properties. The main difference is in the mag- 
nitude of conductance - for parallel DQDs it is roughly 
two times larger than for serial DQDs. This is contrary 
to the strong coupling regime where, for example, in the 
parallel geometry the orbital Kondo phenomenon arises 
due to the interference effects, while for DQDs coupled in 
series the orbital Kondo effect is destroyed^ Finally, we 
also note that in the case of detuned levels, |e2 — £i | 3> 0, 
transport properties of double quantum dots coupled in 
parallel are qualitatively similar to those of single multi- 
level quantum dots^ 
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FIG. 9: (Color online) The differential conductance (a) and 
Fano factor (c) in the parallel magnetic configuration, and 
the TMR (b) as a function of the bias voltage V and the 
position of the dots' levels e = ej = £2 for double quantum 
dots coupled in parallel. The parameters are the same as in 
Fig. [2] with r rj = T/2, for r = L,R, j = 1,2, and T = 0.1 
meV. The transport regime where the Fano factor is divergent 
is marked with a thick black line. 



VI. T-SHAPED DOUBLE QUANTUM DOTS 

An interesting situation occurs when only one of the 
two quantum dots is coupled to external leads, while the 
other one is decoupled. In such T-shaped systems trans- 
port takes place through the molecular states of the dou- 
ble dot, although only the first dot is directly coupled 
to the leads. This may lead to new transport behavior, 
especially to a large super-Poissonian shot noise and an 
enhanced TMR, as we show in the sequel. 

In Fig. [10] we present the current, differential conduc- 
tance, the TMR and the Fano factor as a function of the 
bias voltage for the case when the first dot is coupled 
to the leads while the second dot is decoupled. Spin- 
dependent transport through single quantum dots have 
already been extensively studied. However, the hopping 
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FIG. 10: (Color online) The bias dependence of the current 
and differential conductance (a), the TMR (b), and the Fano 
factor (c) in the parallel (solid line) and antiparallel (dashed 
line) configurations for T-shaped double quantum dots. The 
parameters are: r r i = T/2, and F r 2 = for r = L,R, with 
r = 0.1 meV, and e x = 1 meV, e 2 = meV, k B T = 0.18 
meV, U = 2 meV, U' = 1 meV, t = 0.2 meV, and p = 0.5. 
The dotted curves present the TMR (b) and the Fano factor 
in the parallel configuration (c) calculated for F r 2 = r r i. 



t between the two dots may modify the transport prop- 
erties significantly, giving rise to novel behavior. First 
of all, we note that due to the finite t, transport takes 
place through molecular states of the double dot. This 
in turn leads to additional steps in the i-v curves and 
peaks in the differential conductance, see Fig. fTUT a). as 
compared to tunneling through one single-level quantum 
dots, where for spin-degenerate level the current exhibits 
only two steps. The difference between the currents in the 
parallel and antiparallel configurations gives rise to the 
TMR which is shown in Fig.QUb). The TMR displays an 
oscillatory-like behavior with increasing the bias voltage. 
In addition, at the threshold for sequential tunneling, the 
TMR becomes increased above the Julliere TMR^i This 
result is rather counterintuitive, as the Julliere TMR is 



FIG. 11: (Color online) The bias dependence of the current, 
differential conductance (a), TMR (b) and the Fano factor 
(c) in the parallel (solid line) and antiparallel (dashed line) 
configurations. The parameters are the same as in Fig. [10] 
with £2 = —4 meV. The dotted curves present the TMR (b) 
and the Fano factor in the parallel configuration (c) calculated 

for r r2 = r ra . 



characteristic of a single ferromagnetic tunnel junction 
and one would expect that an additional object (quan- 
tum dot) between the leads should decrease the TMR. 
This was shown for single quantum dots with spin degen- 
erate levels and symmetric couplings to the leads where 
the TMR was found to take at most the Julliere valued 
In principle, the TMR can be enhanced above the Jul- 
liere value when only one spin component of the dot takes 
part in transport. Such spin selection may be achieved 
upon applying an external magnetic field or due to finite 
exchange interaction between spins in the doti 46 i 86 In the 
case of T-shaped DQDs, we find that the enhancement of 
the TMR is associated with increased occupation of the 
second dot in the antiparallel configuration, as compared 
to the parallel. Because the second dot is decoupled from 
the leads, this increases the difference between the cur- 
rent in the two configurations, yielding the TMR larger 
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than the Julliere value. In addition, finite occupation of 
the decoupled dot also leads to large current fluctuations, 
as the rate for tunneling between the second dot is much 
slower than that for tunneling between the first dot and 
the leads. This in turn gives rise to super-Poissonian shot 
noise, as displayed in Fig. riffl e). The super-Poissonian 
shot noise is present in both magnetic configurations. It 
is also slightly larger in the parallel configuration than in 
the antiparallel one, which is due to an additional contri- 
bution to the noise coming from spin-dependent tunnel- 
ing between the dots and the leads^ We note that the 
super-Poissonian shot noise in T-shaped dots has also 
been reported in the case of nonmagnetic leads^i The 
two dots were however rather weakly coupled to each 
other so that no molecular states were formed, contrary 
to the case considered here. 

The enhancement of the TMR above the Julliere value 
and the large super-Poissonian shot noise are present 
approximately at the threshold for sequential tunneling. 
In this transport regime the current is mainly mediated 
through the charge states of the decoupled dot. The 
aforementioned effects should therefore become washed 
out if there were a finite coupling between the second dot 
and the leads. This is shown in Fig. [Tfflb) and (c) where 
we depict the TMR and Fano factor in the parallel con- 
figuration calculated for r r 2 = T r i (see the dotted curves 
in Fig. [TDD. The TMR larger than the Julliere TMR 
accompanied by super-Poissonian shot noise is thus an 
indication that transport takes place through quantum 
dot coupled to another dot, which is not connected to 
the leads directly. 

Figure [Til displays the bias dependence of the current, 
TMR and the Fano factor calculated in the case when 
at equilibrium the second dot is fully occupied while the 
first dot is empty. At low bias voltage, the system is in 
the Coulomb blockade and the double dot is occupied 
with two electrons in the second dot with unit probabil- 
ity. The current flows then only due to elastic cotunneling 
processes. This gives rise to the Julliere TMR and Pois- 
sonian shot noise, see Fig. [TTT b) and (c). With increasing 
the bias voltage, close to the threshold for sequential tun- 
neling, the occupation of the second dot is slightly low- 
ered at the cost of a finite occupation of another states. 
This opens the system for the cotunneling-assisted se- 
quential and sequential tunneling processes. Neverthe- 
less, because the double dot is still mainly occupied by 
two electrons on the decoupled dot, this gives rise to ex- 
tremely large current fluctuations and Fano factors, see 
Fig. [TlTc). This super-Poissonian shot noise can be con- 
siderably reduced once the second dot becomes coupled 
to the leads, see the dotted curve in Fig. [TTT ch 

Finally, in Fig. rT2] we present the differential conduc- 
tance in the parallel configuration and the TMR as a 
function of the bias voltage and the position of the sec- 
ond dot level. The Coulomb blockade regimes are clearly 
visible in Fig. [T27 a). whereas the different behavior of the 
TMR depending on the transport regime is presented in 

Fig. Hub). 
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FIG. 12: (Color online) The differential conductance (a) in 
the parallel magnetic configuration, and the TMR (b) as a 
function of the bias voltage V and the position of the second 
dot level e = £2 for T-shaped double quantum dots. The 
parameters are the same as in Fig. 1101 



VII. CONCLUSIONS 

We have analyzed the spin-dependent transport 
through systems built of two strongly coupled quantum 
dots which are weakly connected to external ferromag- 
netic leads. The considerations were based on the real- 
time diagrammatic technique which allowed us to deter- 
mine transport properties in the sequential and cotunnel- 
ing regimes in a fully systematic way. In particular, we 
have analyzed the current, differential conductance, shot 
noise and the TMR for different geometries of the double 
dot system, including the serial and parallel couplings, as 
well as the T-shaped systems. We have also discussed the 
main differences in transport characteristics correspond- 
ing to different DQDs geometries, which may be helpful 
in determining the geometry of coupled quantum dots in 
experiments. 

In the case of double quantum dots coupled in series, 
we have found an interesting dependence of the TMR 
on the occupation number of the DQD. Furthermore, 
the super-Poissonian shot noise in the Coulomb blockade 
regimes have been observed. On the other hand, when 
the levels of the DQD were detuned, transport character- 
istics revealed the Pauli spin blockade effects. We have 
shown that the leakage current, observed experimentally 
in the blockade regime, 6 results from the interplay of co- 
tunneling processes which flip the spin in the DQD and 
make the sequential tunneling possible. Thus, the current 
in the spin blockade flows due to cotunneling and spin- 
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flip cotunncling-assisted sequential tunneling processes. 
This mechanism is associated with pure nature of tunnel- 
ing processes and is thus relevant for both nonmagnetic 
and ferromagnetic leads. In addition, we have also shown 
that the shot noise in the Pauli spin blockade regime be- 
comes super-Poissonian, while outside the blockade it is 
sub-Poissonian. 

For double quantum dots coupled in parallel, transport 
characteristics were found to be qualitatively similar to 
those in the case of DQDs coupled in series. The main 
difference is in the magnitude of conductance - for par- 
allel DQDs it is roughly two times larger than for serial 
DQDs. 

In addition, we have also analyzed the case when the 
first quantum dot is coupled to the leads while the sec- 
ond one is completely decoupled. In such T-shaped sys- 
tems we have found a large super-Poissonian shot noise at 
the threshold for sequential tunneling and the TMR en- 
hanced above the Julliere value. These effects are associ- 
ated with an increased occupation of the decoupled quan- 
tum dot, and become washed out once there is a finite 
coupling between the second dot and the leads. Thus, the 
enhanced TMR together with large super-Poissonian shot 
noise may be an indication that transport takes place 
through quantum dot which is side-coupled to another 



dot disconnected from the leads. 
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APPENDIX A: DETAILS OF NUMERICAL 
CALCULATIONS 

In order to calculate the transport properties in the 
sequential and cotunncling regimes, it is necessary to 
find the elements of the respective first-order and second- 
order self-energy matrices. This can be done using the 
respective diagrammatic rules i 43 ' 75 ! 76 Here, as an exam- 
pie, we present the contribution coming from W^/ N \ x '/m 
where N is the charge state of the double dot. It is given 

by 



W, 



(2) 



X(N), X '(N) 



-2^EEEE{l^l^ t ^^lx")| 2 |(x"l^ t 4 CT ^lx')| 2 [irt^x" -zx)<r y {s x „ -e x .) 



r,r' j.J' a, a' X " 



Ir'i' ( £ X" ~ £ x') B 2ri( £ x" ~ £ x) ~ B 2^jr<i< i £ X" ~ £ X' £ X' ~ £ x) 



Ir'j' y^X °X ) ± -'2rj\ c -X 



2rjr'j' \°X 



+ \(x\U^dlu\x")\ 2 \{x"\^d ri7 ,U\ X >)\ 2 \^ i {e x -e x „)Bif 3 ,{e x ,-e x n) 

x' 



+rt+(e x , - e x „)B^{e x - e x „) - B££, + (e x - e x „,£ x - e 

+ E(xl^ t 4^lx , '0(x'1^ t ^' CT ^lx')(x|c/ t ^ CT ^lx'')(x''l^ t 4^lx') x 



x'" 



e-y + e. 



B lrjr'f( S X" £ X'-> £ X £ X') B lrjr'f( 6 X £ X'"> £ X £ x') } > ( A ^) 



where 7 ^ ± (e) = r^-/ ± (e - Hr)/{2ir), with /+ being the 
Fermi distribution function, /~ = 1 — / + , fi r denoting 
the electrochemical potential of the lead r, and 

Bgffi (e, e') = J durtf ( U - e<) . 

We note that having found all the first-order and second- 
order self-energy matrices, we are in principle able to 



calculate transport through arbitrary number of differ- 
ent orbital levels coupled to each other and to external 
leads in an arbitrary way. As far as numerical details 
are concerned, for systems consisting of larger number of 
orbital levels, in calculations we make use of the block 
structure of the initial Hamiltonian in the charge space, 
and perform unitary transformation of the Hamiltonian 
and the relevant local operators for each block separately. 
In addition, we also store the respective matrix elements 
in blocks labelled by charge quantum numbers. 
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